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I.  INTRODUCTION  AND  PERSPECTIVE 


It  is  well  known  that  the  inevitable  presence  of  inaccuracies  in  the 
mathematical  models  of  physical  systems  used  for  control  design  can  severely 
compromise  the  resulting  performance.  These  modeling  uncertainties  can  be 
broadly  categorized  as  ( 1 )  parametric  (structured)  uncertainty  resulting  from 
the  variations  in  the  parameters  of  the  low  frequency,  finite  dimensional 
model  and  (2)  high  frequency  unstructured  uncertainty  resulting  from  the 
negligence  of  unmodeled  dynamics.  A  "robust*  control  design  is  that  design 
which  maintains  closed  loop  system  stability  and  performance  even  in  the 
presence  of  these  perturbations. 

The  published  literature  on  this  "robust"  stabilization  and  performance 
problem  can  be  viewed  from  two  perspectives,  namely,  (1)  frequency  domain 
viewpoint,  and  (2)  time  domain  viewpoint.  In  the  frequency  domain  viewpoint, 
the  system  is  assumed  to  be  described  by  transfer  function  representation, 
thereby  making  extensive  use  of  "polynomial"  theory,  whereas  the  time  domain 
treatment  is  naturally  suited  to  systems  represented  by  state  space 
description,  thereby  making  extensive  use  of  "matrix"  theory. 

In  the  transfer  function  based  frequency  domain  representation,  the 
"structured"  uncertainty  is  manifested  in  the  form  of  variations  in  the 
coefficients  of  the  numerator  and  denominator  polynomials,  whereas  the 
"unstructured"  uncertainty  is  reflected  in  terms  of  the  specification  of  a 
norm  bounded  frequency  dependent  function.  On  the  other  hand,  in  the  state 
space,  time  domain  framework,  "structured"  uncertainty  is  manifested  in  the 
form  of  variations  in  the  elements  of  the  system  matrices  while  "unstructured" 
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uncertainty  takes  the  form  of  a  norm  bounded  perturbation  matrix  as  affecting 
the  system. 

In  this  research,  we  employ  both  frequency  domain  viewpoint  as  well  as 
time  domain  viewpoint  in  treating  the  problem  of  stability  robustness,  with 
particular  emphasis  on  "structured"  uncertainty  when  treated  alone.  Towards 
this  direction  we  present  results  on  both  "analysis"  and  on  "design" 
techniques  for  linear  systems  with  structured  uncertainty.  However,  realizing 
that  concentrating  on  only  one  type  of  uncertainty  (either  structured  or 
unstructured)  in  the  problem  formulation  is  inadequate,  since  most  of  the 
applications  of  interest  have  both  types  of  uncertainty  present 
simultaneously,  we  also  treat  the  combined  uncertainty  problem  and  present 
some  interesting  results  on  this  aspect,  again  from  both  time  domain  and 
frequency  domain  viewpoints. 

The  report  is  organized  as  follows:  Section  II  considers  the  analysis 
aspect,  in  which  a  brief  review  of  "polynomial"  and  "matrix"  approaches  in  the 
stability  robustness  analysis  of  linear  systems  is  given.  Then  in  subsections 
2.2  and  2.3,  two  methods  are  presented,  one  for  the  polynomial  case  and  the 
other  for  state  space  case,  to  improve  the  stability  robustness  bounds  under 
dependent  structured  uncertainty.  Section  III  addresses  the  synthesis  aspect 
and  presents  a  time  domain  control  design  method  for  robust  stability  and 
regulation  in  subsection  3.1,  and  in  3.2  a  frequency  domain  based  design 
method  is  given  to  simultaneously  stabilize  a  class  of  non-minimum  phase 
systems.  Section  IV  considers  the  combined  structured  and  unstructured 
uncertainty  problem  and  presents  both  time  domain  based  and  frequency  domain 
based  design  methods  in  subsections  4.1  and  4.2,  respectively.  Finally, 
Section  V  offers  some  concluding  remarks  and  recommendations  for  future 


research. 
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II .  ANALYSIS  FOR  STABILITY  ROBUSTNESS  OF  LINEAR  SYSTEMS 
WITH  STRUCTURED  UNCERTAINTY 

In  this  section,  we  present  some  new  results  on  the  stability  robustness 
analysis  for  linear  systems  with  structured  uncertainty  described  by  both 
transfer  function  representation  and  state  space  representation.  In  the 
former  case,  the  "structured"  uncertainty  is  manifested  in  the  form  of 
variations  in  the  coefficients  of  the  numerator  and  denominator  polynomials  of 
the  transfer  function,  while  in  the  latter  case,  it  is  reflected  in  terms  of 
the  variations  in  the  elements  of  the  system  matrices.  Accordingly,  the 
stability  robustness  analysis  in  transfer  function  formulation  is  addressed 
using  the  Hurwitz  "polynomial"  theory  while  that  of  the  state  space 
formulation  is  addressed  using  "matrix"  theory.  Thus  we  label  the  approaches 
for  the  former  case  as  "polynomial"  approaches  and  those  for  the  latter  case 
as  "matrix"  approaches.  Also,  in  each  of  these  perspectives,  the  robustness 
analysis  problem  can  be  either  in  the  form  of  checking  stability  of  the  family 
of  polynomials  (or  matrices)  given  the  bounds  on  the  parameters  ("Hurwitz 
invariance  tests")  or  in  the  form  of  determining  the  allowable  bounds  on  the 
parameters  to  maintain  stability  ("stability  robustness  measures").  Before 
presenting  the  two  new  results,  one  from  each  perspective,  in  subsection  2.1, 
we  briefly  review  the  available  literature  on  this  aspect. 
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2.1  BRIEF  REVIEW  OF  "POLYNOMIAL"  AND  "MATRIX"  APPROACHES 
Polynomial  Approaches:  In  the  framework  of  continuous  time  uncertain  systems 
represented  by  transfer  function  polynomials,  the  robustness  issue  is  cast  as 
a  problem  of  testing  the  Hurwitz  invariance  of  a  faml 1 v  of  polynomials 
parametrized  by  an  uncertain  parameter  vector  3  6  R1,  which  we  denote  as 
f (s , q) .  The  different  methods  available  for  this  problem  are  dependent  on  the 
characterization  of  the  "structure"  of  the  uncertainty;  i.e.,  the  way  the 
coefficients  of  the  polynomial  depend  on  the  real  parameters  q1 
(i  -  1,2, . . . ,r) .  This  characterization  can  be  divided  into  four  categories: 

Category  PI.  Independent  variations:  In  this  case,  the  coefficients 
vary  independently  of  each  other.  For  example, 

f(s,q)  -  s3  +  3q1s2  +  4q2s  +  5,  Si  <  <U  <  5i»  i  -  1,2. 

The  family  of  polynomials  generated  in  this  category  is  called  an  "interval 
polynomial"  family. 

Category  P2.  Linear  dependent  variations:  For  this  case,  the 
coefficients  vary  linearly  on  the  uncertain  parameters  q1.  For  example, 
f(s,q)  -  s5  +  3q1s*  +  4q1s3  +  2q2sz  +  q2s  +  4q3. 

In  addition,  the  type  of  polynomial 

f(s,g)  -  s3  +  (1  -  q)s2  +  (3  -  q)s  +  (3  -  q) 
is  said  to  be  "affinely  linear  dependent."  The  family  of  polynomials 
generated  in  this  category  is  a  "polytope"  of  polynomials. 

Category  P3.  Multilinear  case:  In  this  case,  the  coefficients  have 
product  terms  in  the  parameters  but  are  such  that  they  are  linear  functions  of 
a  single  parameter  when  others  are  held  constant.  One  example  is 
f(s,g)  -  s3  +  2q1q2s2  +  4q2q3s  +  Sq^q., 
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Category  P4.  Nonlinear  case:  The  coefficients  are  nonlinear  functions 
of  the  parameters.  For  example, 

f(s,g)  -  s3  +  2q?qzs2  +  (4qz  +  3q3)s  +  5 
With  this  classification  in  mind,  we  can  now  present  some  results  available 
for  these  cases.  In  [1],  Guiver  and  Bose  consider  the  "interval"  polynomials 
and  derive  a  maximal  measure  of  robustness  for  quartics  (n  <  4,  where  n  is  the  ! 

degree  of  the  polynomial).  Perhaps  the  best  known  and  most  significant  result 
for  the  "interval"  polynomial  stability  is  the  one  by  Kharitonov  [2] .  This 
result  shows  that  of  the  2n  extreme  polynomials  formed  by  restricting  the  j 

coefficients  to  the  endpoints  of  the  allowable  range  of  variation,  stability 
can  be  determined  by  examining  only  four  special  members  of  this  set 
(independent  of  the  degree  n) .  It  may  be  noted  that  Kharitonov’s  result  is  a 
necessary  and  sufficient  condition  for  the  stability  testing  of  interval 
polynomials  but  becomes  a  sufficient  condition  for  the  case  of  a  polytope  of 
polynomials.  Kharitonov's  result  was  introduced  into  the  western  literature 
by  Barmish  [3].  Later  in  [4],  Bose  recasts  the  Kharitonov  analysis  in  a 
system  theoretic  setting,  and  in  [5],  Anderson,  Jury,  and  Mansour  present  a 
refinement  of  Kharitonov's  result  by  showing  that  fewer  than  four  polynomials 
are  required  for  stability  testing  for  degree  n  <  5. 

Some  additional  results  on  interval  polynomials  and  some  early  results 
on  the  polytope  of  polynomials  can  be  found  in  References  [6]  through  [10], 

The  next  important  result  for  the  "polytope"  of  polynomials  is  given  in 
Bartlett,  Hollot,  and  Lin  [11]  in  which  they  provide  necessary  and  sufficient 
conditions  for  checking  stability.  This  result  is  now  known  as  the  "edge" 
theorem.  It  is  in  a  way  a  Kharitonov- like  result  for  the  polytope  of 
polynomials,  in  the  sense  that,  instead  of  the  vertices  of  a  hypercube,  one 
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has  to  check  the  "edges"  of  the  polytope  to  ascertain  stability.  However, 
this  result  suffers  from  a  "combinatorial  explosion"  problem.  That  is,  the 
number  of  "edges"  to  check  increases  significantly  with  increase  in  the  number 
of  uncertain  parameters.  The  most  recent  result  for  this  problem,  which 
circumvents  somewhat  the  "explosion"  problem,  is  given  by  Barmish  [12].  Here 
a  special  function  called  "robust  stability  testing  function,"  H (6)',  is 
constructed  and  robust  stability  is  assured  if  and  only  if  this  function 
remains  positive  for  all  5  >  0. 

In  subsection  2.2,  we  present  a  result  which  reduces  the 
conservatism  of  the  Kharitonov  test  (as  a  sufficient  condition)  in  the  testing 
of  Hurwitz  invariance  of  a  family  of  polynomials  with  dependent  coefficients. 
This  dependency  can  belong  to  categories  P2,  P3,  or  P4. 

Matrix  Approaches:  In  the  framework  of  continuous  time  uncertain 
systems  described  by  state  space  representation,  the  robustness  issue  is  cast 
as  a  problem  of  testing  the  stability  of  a  family  of  real  n  X  n  matrices, 
parametrized  by  an  uncertain  parameter  vector  cj  e  Rr ,  which  we  denote  as  A(<j)  . 
Specifically,  one  can  write  the  uncertain  system  as 

x(t)  -  A(a)x(t)  -  [A0  +  E(a)]x(t) 

where  A0  is  an  asymptotically  stable  matrix.  The  different  methods  available 
for  this  problem  are  dependent  on  the  characterization  of  the  perturbation 
matrix  E(£) .  If  only  a  norm  on  the  perturbation  matrix  E,  |  |e|  |,  is  known, 
and  bounds  on  it  for  stability  of  A0  +  E  are  derived,  then  it  is  labeled  as 
the  "unstructured"  uncertainty  problem.  Some  of  the  results  for  this  case  are 
reported  in  References  [13]  through  [21],  However,  in  this  research  we 
concentrate  mainly  on  the  "structured"  uncertainty  case.  For  this  case,  as  in 


6 


the  polynomial  setting,  the  different  methods  available  are  again  influenced 
by  the  characterization  of  the  "structure"  of  the  uncertainty;  namely,  the  way 

the  elements  of  E  depend  on  the  real  parameters  qA  (i  -  1,2 . r) .  Again,  we 

can  consider  the  following  categories: 

Category  Ml.  Independent  variations:  In  this  case  the  elements  of  E 
vary  independently,  i.e., 

-ij  -  Eij  -  Eu  ^or  *-•■!  “  1*2 . n 

Another  way  of  representing  this  situation  is 
r-n2 

£(s>  -  ^iEi 

i-i 

where  E1  are  known  constant  matrices  with  only  one  nonzero  element,  at  a 
different  location  for  different  i. 

This  family  of  matrices  is  the  so-called  "interval  matrix"  family. 

Cfltegpcy  M2. _ lingflt  dependent  variation:  For  this  case,  the  elements 

of  E  vary  linearly  with  the  parameters  qt.  Thus,  we  can  write 
r 

E(q)  -  l  qiEi 
i-1 

where  Et  are  constant,  known  real  matrices.  This  family  is  termed  the 
"polytope"  of  matrices. 

Cflt-EgQry.ML _ General  variation:  For  this  case 


E(a>  - 


r 


I 


i-1 


q^Ca) 


Unlike  the  polynomial  case,  the  matrix  case  is  decidedly  more  difficult, 
as  evidenced  by  the  lack  of  strong  results.  In  other  words,  even  though  there 
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are  "derivable"  necessary  and  sufficient  conditions  available  [22]  for  any 
general  type  of  dependency,  these  conditions  are  not  really  that  helpful  from 
tne  practical  application  point  of  view  because  the  conditions  are  not  simple 
or  computationally  tractable.  It  can  be  safely  said  that  at  the  present  time 
there  are  no  simple,  easily  and  finitely  computable  necessary  and  sufficient 
conditions  available  even  for  the  simplest  "interval"  matrix  case.  For  some 
special  cases  of  interval  matrices,  the  results  of  References  [23] -[26]  are 
helpful.  One  natural  line  of  thought  for  solving  this  problem  (say  "interval" 
matrix  stability  problem)  would  be  to  transform  the  interval  matrix  stability 
problem  to  the  "interval"  or  "polytope"  of  polynomial  problem  and  then  apply 
the  results  discussed  previously.  However,  extreme  caution  is  warranted  in 
this  approach  in  view  of  the  abundance  of  counter-examples,  as  discussed  in 
References  [27] -[29].  Also,  this  approach  is  necessarily  conservative  (even 
when  it  works)  as  shown  by  Yedavalli  in  Reference  [30],  In  view  of  these 
difficulties,  much  of  the  literature  on  the  matrix  family  stability  problem  is 
dominated  by  results  which  provide  "sufficient"  conditions  for  stability. 

In  the  framework  of  testing  the  stability  of  "interval"  matrices,  Heinen 
[31]  and  Argoun  [32]  present  results  which  use  Gershgorin's  diagonal  dominance 
concept.  Xu  [33]  imposes  certain  leading  principal  minor  conditions  and  Shi 
and  Gao  [34]  assume  "symmetry"  (thereby  solving  a  special  case).  In  the 
framework  of  obtaining  measures  of  robustness  (the  results  of  which  can  also 
be  applied  to  test  stability),  Patel  and  Toda  [16],  Yedavalli  [35],  and 
Barmish  and  DeMarco  [36]  use  the  Lyapunov  method.  Presently,  Juang  et.  al. 
[37]  and  Qiu  and  Davison  [38]  used  frequency  domain  arguments  to  present 
sufficient  bounds  for  stability  of  the  interval  matrix  problem. 


8 


For  the  case  of  the  "polytope”  of  matrices,  Zhou  and  Khargonekar  [39] 
extend  the  Lyapunov  method  used  by  Yedavalli  [19] -[20]  to  give  sufficient 
conditions  for  stability  and  show  that  the  independent  variation  case  of 
Yedavalli  becomes  a  special  case  of  their  result.  Thus  they  get  better  bounds 
for  the  linear  dependent  variation  case.  In  the  present  research,  in  subsection 
2.3,  we  present  new  results  on  obtaining  sufficient  conditions  for  the  linear 
dependent  variation  case  as  well  as  for  the  nonlinear  case,  using  the 
frequency  domain  results  of  Reference  [37]. 

2.2  REDUCED  CONSERVATISM  IN  THE  TESTING  OF  HURWITZ  INVARIANCE 
OF  UNCERTAIN  POLYNOMIALS  WITH  DEPENDENT  COEFFICIENTS 

In  the  classical  stability  problem  of  linear  time- invariant  systems,  one 
is  often  required  to  decide  whether  the  zeros  of  a  given  polynomial  f(s)  lie 
in  the  strict  left  half  plane.  The  renowned  Routh-Hurwitz  criterion  [40]  has 
proven  to  be  a  valuable  tool  for  this  purpose.  More  recently,  however,  much 
attention  has  been  devoted  to  so-called  robustness  analysis.  Given  the 
motivation  that  system  parameters  are  often  uncertain  or  changing  with  time, 
the  fixed  polynomial  f(s)  is  replaced  by  a  family  of  polynomials 
(f(s,q):  q  E  Qj  where  Q  is  a  compact  set  and  one  is  faced  with  the  problem  of 
determining  whether  f(s,q)  is  stable  for  all  q  €  Q.  In  view  of  the  fact  that 
applying  the  Routh-Hurwitz  criterion  to  each  member  of  the  set  Q  is 
impractical,  a  number  of  alternative  methods  are  proposed  as  discussed  in  the 
previous  section.  Perhaps  most  notable  are  the  results  of  Kharitonov  [2], 
which  apply  to  the  special  case  obtained  when  the  coefficients  of  f(s,q)  vary 
independently.  In  this  situation,  f(s,q)  is  Hurwltz  invariant  over  Q  if  and 
only  if  four  specially  constructed  "bounding  polynomials"  are  stable. 
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However,  in  many  control  problems,  the  coefficients  of  resulting  uncertain 
polynomials  f(s,q)  usually  depend  on  each  other  and  the  Kharitonov  test 
becomes  a  sufficient  condition;  as  such,  it  becomes  conservative. 

In  this  note  we  propose  a  new  criterion  to  reduce  the  conservatism  of 
the  Kharitonov  test  for  Hurwitz  invariance  of  uncertain  polynomials.  It  can 
be  viewed  as  a  generalized  version  of  Kharitonov's  results  and  it  is  typically 
suitable  for  dealing  with  uncertain  polynomials  having  dependent  coefficients. 
This  result  is  organized  as  follows:  In  subsection  2.2.1,  we  formally  define 
uncertain  polynomials  and  provide  some  lemmas  to  be  used  in  later  sections.  In 
subsection  2.2.2,  we  state  the  main  results.  Some  illustrative  examples  are 
given  in  subsection  2.2.3. 

2.2,1  Preliminary  Definitions  and  Lemmas 

Definition  2.1:  An  n-th  order  uncertain  polynomial  is  defined  by 
f(s,q)  -  aQ(q)sn  +  o^qjs"'1  +  ...  +  ajq)  (2.1) 

where  q  is  the  uncertainty  index  vector,  Q  is  a  compact  subset  of  Rp,  which  is 
the  index  set  for  q,  at:  Q  -»  R  are  real  continuous  coefficient  functions. 
f(s,q)  is  said  to  be  monic  if  a0(q)  ■  1.  f(s,q)  is  said  to  have  independent 
uncertain  coefficients  if  q  -  (q0  q:  •  •  •  qn)  and  aA(q)  -  qA  for  all 
i  -  0,1,..., n.  f(s,q)  is  said  to  have  dependent  uncertain  coefficients  if 
there  exist  at(q)  -  ^(oj(q),q)  for  i  /  J. 

Definition  2.2:  An  uncertain  polynomial  f(s,q)  as  in  (2.1)  is  said  to 
be  Hurwitz  invariant  over  Q  if  for  each  q  €  Q  all  the  zeros  of  f(s,q)  lie  in 
the  strict  left  half  plane. 

The  following  lemma  is  an  immediate  consequence  of  the  Routh-Hurwitz 
theorem.  Therefore,  no  proof  is  provided. 
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Lamina  2.1:  An  uncertain  polynomial  f(s,q)  as  in  (2.1)  with  a0(q)  >  0  is 


Hurwitz  invariant  if  and 

only  if  the 

Hurwitz 

testing 

matrix 

“i(q) 

<*3(q) 

o3(q> 

0 

0 

o0<q) 

«2<q) 

«*<q) 

0 

0 

H(q)  - 

0 

“i(q) 

<*3(q) 

a5(q> 

0 

0 

«0<q) 

“2(q) 

«**(q) 

0 

0 

0 

«i(q) 

a3(q) 

<*s(q> 

(illustrated  above  for  n  -  5)  has  positive  leading  principal  minors  for  all 

q  6  Qn- 

Lemma  2.1  provides  a  solid  theoretical  basis  for  testing  Hurwitz 
invariance  of  uncertain  polynomials.  However,  when  an  uncertain  polynomial 
has  high  order  and  its  index  set  Q  contains  an  infinite  number  of  elements,  it 
is  impractical  to  apply  the  lemma.  A  computationally  feasible  test  for 
Hurwitz  invariance  is  needed. 

Definition  2.3:  Four  bounding  polynomials  of  an  uncertain  polynomial 


f(s,q)  as  in  (2.1)  are  defined  as  follows: 

fx(s)  ■  a0sn  +  ajs"'1  +  b2sn~2  +  b3s“'3  +  a4sn'*  +  ...  (2.2) 

f2(s)  ■  b0sn  +  bjs"'1  +  a2sn‘2  +  a3s”*S  +  b4s"'*  +  ...  (2.3) 

f3(s)  ■  a0sn  +  bjSn"1  +  b2sn*2  +  a3sn'3  +  a4sn*4  +  ...  (2.4) 

f4(s)  ■  b0sn  +  SjS11'1  +  a2sn'2  +  b3s“'3  +  b4s”*4  +  ...  (2.5) 

where  for  all  i  -  0,1,..., n, 

ax  -  min{a1(q):  q  e  Q) ,  and 
bt  -  max(o1(q) :  q  6  Q) . 
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Lemma  2 . 2  (the  Kharitonov  theorem;  see  Reference  [2]  for  proof):  An 
uncertain  polynomial  f(s,q)  as  in  (2.1)  is  Hurwitz  invariant  over  Q  if  its 
four  bounding  polynomials  are  Hurwitz.  If,  in  addition,  all  the  coefficients 
of  f (s,q)  are  independent,  then  f(s,q)  is  Hurwitz  invariant  over  Q  if  and  only 
if  its  four  bounding  polynomials  are  Hurwitz. 

Remark:  In  Reference  [2],  only  monic  uncertain  polynomials  are 
considered.  However,  with  minor  modification,  the  Kharitonov  theorem  can  be 
stated  in  a  form  as  in  Lemma  2.2  when  a0(q)  is  not  identity. 

Note  that  in  most  control  problems  we  are  actually  dealing  with 
uncertain  polynomials  having  dependent  coefficients,  and  it  is  often  too 
conservative  for  Hurwitz  invariance  to  require  four  Hurwitz  bounding 
polynomials.  Our  main  goal  in  this  section  is  to  provide  a  methodology  to 
reduce  this  conservatism. 

Mfrlh  Results 

Theorem  2.1  (see  Appendix  A  for  proof):  An  uncertain  polynomial  f(s,q) 
as  in  (2.1)  is  Hurwitz  invariant  over  Q  if  and  only  if  there  exist  positive 
continuous  functions  kj(q)  >  0,  k2(q)  >  0  for  all  q  e  Q  such  that 

f’(s,q)  -  a’(q)an  +  a[(q)a^1  +  ...  +  c*;<q),  qeQ  (2.6) 

where 

(ki(q)«i(q).  i  even; 

MqJVq).  1  odd- 

is  Hurwitz  invariant  over  Q. 

Theorem  2.1  provides  a  way  to  transform  an  uncertain  polynomial  to 
another  without  changing  its  Hurwitz  invariance.  Hence,  combining  Theorem  2.1 
and  Lemma  2.2,  we  immediately  have  the  following  theorem. 
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Theorem  2.2:  When  positive  continuous  functions  kx(q)  >  0,  k^q)  >  0 
for  all  q  e  Q  transform  an  uncertain  polynomial  f(s,q)  to  an  uncertain 
polynomial  f'(s,q)  as  in  (2.6),  then  f(s,q)  is  Hurvitz  invariant  over  Q  if  the 
four  bounding  polynomials  of  f'(s,q)  are  Hurwitz.  If,  in  addition,  all  the 
coefficients  of  f'(s,q)  become  independent,  then  f(s,q)  is  Hurwitz  invariant 
over  Q  if  and  only  if  the  four  bounding  polynomials  of  f'(s,q)  are  Hurwitz. 

Remark :  Theorem  2.2  can  be  viewed  as  a  generalized  version  of  the 
Kharitonov  theorem.  It  provides  some  flexibility  in  applying  the  "four 
bounding  polynomial"  criterion.  By  properly  choosing  kx(q)  and  k2(q) ,  one  can 
reduce  the  conservatism  of  the  Kharitonov  test  for  uncertain  polynomials 
having  dependent  coefficients.  In  the  best  case,  when  kx(q)  and  kj(q)  are 
found  such  that  f'(s,q)  has  independent  coefficients,  the  conservatism  is 
completely  reduced.  If,  however,  such  k2(q)  and  k2(q)  cannot  be  found,  it  may 
still  be  possible  to  partially  reduce  the  conservatism  by  choosing  kx(q)  and 
k2(q)  either  to  relax  the  dependence  of  coefficients  or  to  uniformly  reduce 
the  varying  bounds  of  uncertain  coefficients.  In  the  latter  case,  we  have  the 
following  theorem. 

Theorem  2 . 3  (see  Appendix  B  for  proof):  If  Hurwitz  invariance  of  f(s,q) 
can  be  determined  by  the  Kharitonov  theorem,  it  can  also  be  determined  by 
Theorem  2.2,  provided  kx(q)  and  kj(q)  are  chosen  so  that 

ai  <  al ■  M  <  bi 

where 

a4  -  min{ai(q) :  q  e  Q} , 

bA  -  max(at(q) :  q  €  Q) , 
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-{ 


min{k1(q)a1(q) :  q  6  Q),  i  even; 
min{k2(q)a1(q)  :  q  e  Q) ,  1  odd; 


and 


-I 


max(k1(q)a1(q) :  q  €  Q) ,  i  even; 
oax{k2(q)o1(q) :  q  6  Q),  i  odd. 

Remark  2.3:  Theorem  2.3  indicates  a  way  to  choose  kx(q)  and  k^q)  which 
will  definitely  reduce  the  conservatism  of  the  Kharitonov  test.  If  the 
Kharitonov  test  works,  our  method  also  works.  If  the  Kharitonov  test  fails, 
our  method  may  still  work  because  the  varying  bounds  are  reduced.  Generally 
speaking,  one  should  make  his  own  choice  of  k^q)  and  k2(q)  for  each 
individual  problem.  It  may  also  be  noted  that  the  reduction  in  conservatism 
obtained  by  the  proposed  method  is  in  a  way  similar  in  concept  to  the 
reduction  in  conservatism  obtained  by  state  transformation  by  Yedavalli  and 
Liang  in  References  [41]  and  [35]  In  the  context  of  stability  conditions  for 
interval  matrices. 

Illustrative  Examples 

Example  2.1:  Consider  a  linear  time  invariant  uncertain  system 


x(t)  -  A(q)x(t)  - 


0  3  1 

-1  q  0 

-1  -4  -2 


x(t) ,  q  6  Q  -  [-1,1] 


Determine  whether  the  uncertain  system  is  stable  for  all  q  €  Q. 

This  problem  is  equivalent  to  testing  Hurwitz  invariance  of  the 
following  uncertain  polynomial 

f(s,q)  -  det(sl  -  A(q))  -  s3  +  (2  -  q)s2  +  (4  -  2q)s  +  (2  -  q)  . 
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By  using  the  Kharitonov  theorem,  we  teat  if  the  following  four  bounding 
polynomials  are  Hurwitz: 

f1(s,q)  -  s3  +  s2  +  6s  +  3 
f2(s,q)  -  s3  +  3sz  +  2s  +  1 
f3(s,q)  -  s3  +  3s2  +  6s  +  1 
f4(s,q)  -  s3  +  s2  +  2s  +  3. 

Obviously,  f4(s,q)  is  not  Hurwitz,  hence,  the  Kharitonov  test  fails  in  this 
problem. 

Now,  applying  Theorem  2.3,  we  take  k^q)  -  1  and  kjCq)  -  1/(2  -  q)  and 

obtain 

f'(s.q)  -  s3  +  s2  +  (4  -  2q)s  +  1 

which  has  independent  coefficients.  Since  the  bounding  polynomials  of  f'(s,q) 
f [ (s , q)  -  s3  +  s2  +  6s  +  1 

and 

f2(s,q)  -  s3  +  s2  +  2s  +  1, 

are  Hurwitz,  we  conclude  that  f(s,q)  is  Hurwitz  Invariant  and  the  system  is 
stable  for  all  q  €  Q. 

Example  2.2:  Given  an  uncertain  polynomial 
f(s,q)  -  s4  +  s3  +  2qs2  +  s  +  q,  q  e  [1.5,4], 
determine  if  f(s,q)  is  Hurwitz  invariant. 

It  is  easy  to  see  that  one  of  the  bounding  polynomials, 
f(s)  -  s4  +  s3  +  3s2  +  s  +  4, 

is  not  Hurwitz ,  hence ,  the  Kharitonov  test  also  fails  in  this  example .  Now 
letting  kj(q)  -  1/q  and  kj(q)  -  1,  we  obtain 
f '  (s,q)  -  (l/q)s4  +  s3  +  2s2  +  s  +  1, 
which  has  independent  coefficients.  Since  the  bounding  polynomials 
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fj(s,q)  -  (2/3)s*  +  s3  +  2s2  +  s  +  1 


and 

f2(s,q)  -  (l/4)s*  +  s3  +  2s2  +  8  +  1 
are  Hurwltz,  f(s,q)  is  also  Hurwltz  Invariant. 

Example  2.3:  Given  an  uncertain  polynomial 

f(s,q)  -  s3  +  2qjS2  +  (q2  +  q2)s  +  q^,  qx  6  [1,4],  q2  €  [2,4] 
determine  if  f(s,q)  is  Hurwltz  invariant. 

Since  one  of  the  bounding  polynomials 
f(s)  -  s3  +  2s2  +  3s  +  16 

is  not  Hurwltz,  the  Kharitonov  test  fails.  We  now  choose  kx(q)  -  1  and 
k2(q)  -  l/q:  and  obtain 

f'(s,q)  -  s3  +  2s2  +  (qi  +  q2)s  +  q2,  q2  G  [1,4],  q2  G  [2,4]. 

This  polynomial  still  has  dependent  coefficients.  However,  the  varying  bounds 
of  coefficients  are  uniformly  reduced.  The  four  bounding  polynomials  of 
f' (s.q) , 

f 2 ( ® , q)  -  s3  +  2s2  +  8s  +  4, 

f2(s,q)  -  s3  +  2s2  +  3s  +  2, 

fj(s,q)  -  s3  +  2s2  +  8s  +  2, 

fi(s,q)  -  s3  +  2s2  +  3s  +  4 

are  Hurwltz.  Hence,  f(s,q)  is  Hurwltz  invariant. 

2.3  STABILITY  ROBUSTNESS  MEASURES  FOR  STATE  SPACE  MODELS 
WITH  DEPENDENT  UNCERTAINTY 

In  this  section,  we  move  on  to  the  case  of  the  time  domain  state  space 
representation  and  give  improved  bounds  for  stability  robustness  for  systems 
with  dependent  uncertainty. 
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As  noted  in  section  2.1,  the  aspect  of  developing  explicit  upper  bounds 
on  the  perturbation  of  linear  state  space  systems  to  maintain  stability  has 
received  much  attention,  in  particular,  for  the  problem  of  "structured 
uncertainty."  In  [36],  the  perturbation  matrix  elements  are  taken  to  vary 
independently  of  each  other,  whereas  [39]  treats  the  case  in  which  the 
uncertain  parameters  enter  linearly  into  the  perturbation  matrix.  These  two 
techniques  use  the  Lyapunov  approach  and  can  also  accommodate  time -varying 
(bounded)  perturbations.  On  the  other  hand,  the  techniques  in  [37] -[38]  also 
develop  bounds,  for  independent  variations,  using  a  root  locus  type  argument, 
and  obtain  a  frequency  domain  based  formula  for  the  calculation  of  the  bounds. 
These  techniques  are  basically  aimed  at  time- invariant  perturbations. 

In  this  section  we  present  new  bounds  for  the  dependent  variation  case, 
in  particular  for  the  following  two  cases:  (1)  linear  dependency,  and  (2) 
quadratic  variation  in  a  scalar  uncertain  parameter.  The  proposed  technique 
includes  the  bounds  of  [37] -[38]  as  a  special  case.  The  method  is  illustrated 
with  the  help  of  several  examples. 

Consider  the  state  space,  perturbed  model 

x(t)  -  (A0  +  E)x(t)  (2.7) 

where  A0  Is  an  n  x  n  asymptotically  stable  matrix  and  E  is  a  "perturbation" 
matrix. 

Independent  variations  (Category  Ml):  In  this  case,  the  elements  of  the 
matrix  E  are  assumed  to  be  independent  of  each  other  and  are  bounded  such  that 

Eu  <  lEijL«  “  £ij  and  £  ^  Ma*  cij  (2-8> 

t.J 

Denoting  A  as  the  matrix  formed  with  c1J 

A  -  Kj]  (2.9a) 
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we  write 


A  -  cU.  (2.9b) 

where 

0  <  U#iJ  <  1  (2.9c) 

For  this  case,  the  bounds  on  eiJ  are  obtained  in  [20]  and  [36]  by  a  Lyapunov- 
based  approach  and  in  [37] -[38]  by  a  frequency  domain  based  approach.  In  this 
section,  we  are  interested  in  the  frequency  domain  based  approach,  and  hence 
reproduce  the  expressions  given  in  [37] -[38]  for  the  above  notation.  The 
stability  robustness  bounds  on  e1j  are  given  by  [37] -[38]  ([38]  considers  a 
more  general  structure  also) : 


_  1 _ 

CiJ  <  sup  p[  I  (jwl  -  A0)'1  |UJ 


U.ij  or  £  <  f*lBd 


(2.10) 


w>0 

where  |(*) I  denotes  the  absolute  matrix  (i.e.,  matrix  with  absolute  values  of 
the  elements)  and  p(-)  denotes  the  spectral  radius  of  the  matrix  (•)  and  /ilnd 
denotes  the  bound  for  independent  variations. 

Linear  dependent  uncertainty  (Category  M2):  In  what  follows,  we  consider  the 
case  where  the  uncertain  parameters  in  E  are  assumed  to  enter  linearly,  i.e., 
r 

E  -  l  P, Et  (2.11) 

i-1 


where  Et  are  constant,  prescribed  matrices  and  Pt  are  uncertain  parameters. 

Our  intention  is  to  give  a  bound  on  1^1. 

Ue  now  present  a  bound  on  | Pi |  and  show  that  the  resulting  bound 
specializes  to  (2.10)  for  the  independent  variation  case.  The  proposed  bound 
is  less  conservative  than  (2.10)  when  applied  to  the  situation  in  which  E  is 
given  by  (2.11)  and  yields  the  same  bound  as  in  (2.10)  when  applied  to  the 
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independent  variation  case  (when  E  is  given  by  (2.8)).  This  is  exactly  the 
type  of  situation  that  arises  in  Zhou  and  Khargonekar  [39],  where  they 
consider  the  linear  dependency  case  and  specialize  it  to  the  independent 
variation  case  of  Yedavalli,  [36]. 

Remark:  It  should  be  mentioned  at  the  outset  that  it  is  very  important  to 
distinguish  between  the  independent  variation  case  and  the  dependent  variation 
case  at  the  problem  formulation  stage.  In  the  independent  variation  case  one 
gives  bounds  on  ciJ  (and  consequently  on  e),  whereas  in  the  dependent 
variation  case,  one  gives  bounds  on  1^1.  This  is  particularly  crucial  in 
the  comparison  of  different  techniques.  Proper  comparison  is  possible  only 
when  the  basis,  namely  whether  one  Is  considering  the  dependent  case  or  the 
independent  case,  is  established  beforehand.  It  may  be  noted  that  the 
techniques  aimed  at  the  independent  variation  case  can  accommodate  the 
dependent  variation  situation,  but  at  the  expense  of  some  conservatism; 
whereas  the  technique  aimed  at  the  dependent  case,  while  it  gives  a  less 
conservative  bound  for  that  case,  cannot  accommodate  the  independent  variation 
case  (unless  it  is  established,  additionally,  to  handle  that  situation,  as  is 
done  in  [39]  and  also  later  in  this  section  for  the  new  proposed  bound). 
Theorem  2.4:  Consider  the  system  (2.7)  with  E  as  in  (2.11).  Then  (2.7)  is 
stable  if 

|j®i  I  <  - - - - - for  r  >  1  (2.12a) 

sup  pi  l  |(jwl  -  A0)'1E1|  ] 

w>0  1  i-1  J 

and 


_ i _ 

sup  p[(J«I  -  A0)‘1E1] 
w>0 


jid  for  r  -  1 


(2.12b) 
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where  ni  denotes  the  bound  for  "dependent"  variation  case. 
Proof:  Given  in  Appendix  C. 


It  can  be  shown  that  the  bound  (2.10)  becomes  a  special  case  (2.12a) 
when  one  notes  that  in  the  independent  veriation  case  each  Et  will  contain  a 
single  element  and  is  given  by 


"Hn-D+J  ”  U»ijelej 


(2.13) 


where  e1  is  an  n-dimensional  column  matrix  with  1  in  the  i-th  entry  and  0 
elsewhere.  Note  that  U#1j  is  a  scalar  and  eiej  is  a  matrix. 

Remark :  The  bound  of  (2.12a),  when  specialized  to  the  independent  variation 
case  (i.e.,  when  each  Ei  contains  a  single  element,  at  a  different  place  for 
different  i)  ,  will  be  denoted  by  #ilnd.  Thus,  -  fij  -  fiQ. 

Example  2.4:  Consider 


Ao  " 


-3 

1 


-2 

0 


and  let 


E  -  ft 


1 

1 


1 

1 


(i.e.,  dependent  case) 


Yedavalli[20] 

Hq  [37] 

[38] 

Md  (proposed) 

Zhou,  et.al. [39] 

0.236 

0.329 

0.329 

1.0 

1.0 

If,  instead,  all  the  elements  in  E  are  assumed  to  vary  independently,  then  we 
use 


Ei  “ 

1  0 

0  0 

»  “ 

0  1 

0  0 

.  e3  - 

0  0 

1  0 

.  e4  - 

0  0 

0  1 

- 

. 

. 

.  m 

in  the  expression  (2.12)  and  get 
I  Ei j  I  <  Min<i  "  0  •  329 

which  is,  of  course,  the  same  bound  as  n3  and 

It  may  be  noted  that  the  bound  /id  -  1.0  is  also  a  necessary  and 
sufficient  bound  if  p  is  varied  in  the  positive  direction.  However,  there  was 
no  attempt  mat  a  in  (2.12)  to  incorporate  directional  information. 

Example  2.5:  Consider  the  same  A„  as  in  Example  2.4  and  let 


E 


1 

0 


Yedavalli[20] 

MQ  [38] 

Mj  [37] 

Md  (proposed) 

Zhou,  et.al.[39] 

1.0 

1.52 

1.52 

o 

CM 

2.0 

Example  2.6:  Let  us  consider  the  example  given  in  [39]  in  which  the  perturbed 
system  matrix  is  given  by 


-2  +  kj  0  -1  +  k3 


(A0  +  BKC)  - 


0  -3  +  kj  0 


-1  +  kj  -1  +  kj  -4  +  kj 


Taking  the  nominally  stable  matrix  to  be 


-2 

0 

1 

Ao  “ 

0 

-3 

0 

-1 

-1 

4 

the  error  matrix  with 

kx  and  k2 

as 

the  uncertain  parameters 

Is  given  by 

E  -  k^ 

+  k^ 

where 

1 

0 

1 

0 

0  0 

Ei  “ 

0 

0 

0 

and  E2  " 

0 

1  0 

1 

0 

1 

0 

1  0 

The  following  are  the  bounds  on  |kj|  and  |kj|  obtained  by  [25]  and  the 
proposed  method  (corresponding  to  the  bound  given  by  expression  (2.9)  in 
[25]): 


Yedavalli[20] 

M„  08] 

Mj  137] 

(proposed) 

Zhou,  et.al.[39] 

0.815 

0.875 

0.875 

1.75 

1.55 

It  may  be  noted  that,  in  this  result,  we  do  not  have  counterparts  to  the 
bounds  given  by  expressions  (2.8)  and  (2.7)  in  Zhou  and  Khargonekar ' s  paper. 
Quadratic  Variation  in  a  Scalar  Parameter:  We  now  consider  the  case  in  which 
a  scalar  uncertain  parameter  enters  into  E  in  a  nonlinear  manner,  in 
particular,  as  a  square  term;  i.e.,  we  assume 

EtJ(q)  -  k^q2  (2.14) 

where  q  is  che  uncertain  parameter. 
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This  type  of  structured  uncertainty  model  occurs  in  many  applications; 
for  example,  in  [42]  where  a  large-scale  interconnected  system  has  the  form 


xi  - 


0  1 
-1  -1 


*1  + 


*2  “ 


0  1 
•  1  -1 


x2  + 


which  can  be  written  as 


Ao  “ 


0 

■1 

0 

0 


1 

■1 

0 

0 


0  0 

q2  o 


xi  + 


0  0 
0  -2q 


xl  + 


0 

0 

0 

-1 


0 

0 

1 

-1 


0  0 
0  2q 


0  0 

q2  0 


and  B(q)  - 


0  0 

q2  0 

0  0 

0  -2q 


0 

0 

0 


0 

2q 

0 

0 


(2.15a) 


(2.15b) 


For  a  situation  of  this  type,  the  following  iterative  method  is  proposed  to 
obtain  an  improved  bound  on  | q  | . 


Proposed  Iterative  Method: 

Iteration  1 

In  this  iteration,  we  ignore  the  functional  dependence  and  assume  the 
entries  in  the  perturbation  matrix  vary  Independently.  Accordingly,  we  let 
UslJ  -  1  for  those  entries  in  which  perturbation  is  present  and  zero  for  the 
other  entries.  Then  compute  the  bound  ^  using  the  expression  (2.10).  Let 
the  upper  bound  matrix  be  denoted  by 


K 

i 


- 


(2.16) 
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Knowing  the  elements  of  A,  and  the  corresponding  functional  relationship  of 

1 

the  perturbation  matrix  elements  on  q  in  the  matrix  [E(q)]a  ((•)  denotes  the 

modulus  matrix),  solve  for  the  possible  different  values  of  |q|  and  select 

the  minimum  value  of  |q|.  Let  this  value  of  |q|  be  denoted  by  qa  .  With 

1 

this  value  of  qa  ,  compute  the  matrix  [E(qa  )  ]a  utilizing  the  functional 
1  l 

relationship  in  the  matrix  [E(q)]a.  Then  write 

[E(q.  )]»  a  c2u,2  <2.i7) 

i 

where  c2  is  the  maximum  modulus  element  in  [E(qa  )]a. 

Iteration  2 

Uith  c 2  as  the  left-hand  side  of  Equation  (2.10),  compute  the  bound  jj2 
using  Equation  (2.10)  (with  U#2  replacing  the  matrix  Us) . 

If 

c2  <  (2.18) 

form 

A.  "  ^U.2  (2-19) 

2 

and  go  through  the  exercise  outlined  after  Equation  (2.16)  in  iteration  1  to 

obtain  qB^  ,  which  will  be  greater  than  ra  .  If  e2  {  n2,  the  ra  becomes  the 
2  11 

acceptable  bound. 

Termination:  Repeat  the  iterations  until  no  improvement  in  the  bound  q,, 

l 

(i  -  1,2,...)  is  observed  (say  at  iteration  N)  and  take  q,,,  as  the  acceptable 
bound . 

F.xamnle  2.7:  Consider  the  system  given  in  (2.15). 

Iteration  1 
Let 
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1 


0  0  0 
0  0  1 
0  0  0 
110 


and  compute  ^  -  0.4841.  Form 


0 


0  0 


0 


0.4841  0 
0  0 


0  0.4841 
0  0 


0.4841 


0.4841 


Knowing 


0  0  0  0 

q2  0  0  2q 

0  0  0  0 

0  2q  q2  0 


we  can  solve  for  |q|  as  1)  |q|  -  0.2421;  11)  |q|  -  0.6957.  We  take 


q,,  -  Min[0. 2421,0.6957]  -  0.2421. 

l 

We  then  form  the  matrix  [ECq,  )]„,  l.e., 
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0 

0 

0 

0 

[E(q«)]«  - 

(0.2421)2 

0 

2(0.2421) 

i 

0 

0 

0 

0 

0 

2(0.2421) 

(0.2421)2 

0 

-  0.4841 


0 

0.121 

0 

0 


0 

0 

0 

1 


0 

0 

0 

0.121 


0 

1 

0 

0 


-  caU. 


(Thus  c2  -  0.4841.) 


Iteration  2:  Compute  the  bound  fi2  with  U0  as  the  "structured"  matrix 

2 

and  obtain 

M2  -  0.8921. 

Noting  that  n2  >  c2  (and  that  n2  >  M*) »  we  proceed  further  and  form 


"  M2U,  - 

2  2 


0 

0.1079 

0 

0 


0 

0 

0 

0.8921 


0 

0 

0 


0 

0.8921 

0 


0.1079 


Comparing  A,,,  with  (E(q)]B,  we  can  solve  for  |q|  as  i)  |q|  -  0.3285  (from  qz 
2 

term),  ii)  |q|  -  0.446  (from  2q  term).  We  take 

q,,  -  Min[0. 3285,0.446]  -  0.3285. 

2 

Thus  one  acceptable  range  of  q  is 
0  <  | q  |  <  0.3285. 

One  can  carry  out  these  iterations  further  to  obtain  an  improved  bound  on 
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III.  SYNTHESIS  OF  ROBUST  CONTROLLERS  FOR  LINEAR  SYSTEMS 
WITH  STRUCTURED  UNCERTAINTY 

In  the  previous  section,  the  results  obtained  were  essentially 
"analysis"  tools.  No  attempt  was  made  to  synthesize  controllers.  Of  course, 
an  exercise  in  that  type  of  analysis  is  helpful  before  one  can  design 
controllers  for  a  given  objective.  Aided  with  this  type  of  analysis,  in  this 
section  we  emphasize  the  aspect  of  design  and  present  two  design  methods,  one 
in  the  state  space  framework  and  the  other  in  the  transfer  function  framework. 
In  the  next  section,  the  system  is  described  by  a  state  space  representation 
and  the  aim  is  to  design  a  controller  that  maximizes  the  stability  robustness 
bounds  for  a  given  structured  uncertainty  with  adequate  control  effort 
considerations.  The  design  is  applicable  to  systems  with  time  varying 
perturbations.  In  section  3.2,  the  uncertain  system  is  assumed  to  be 
represented  in  transfer  function  form,  parametrized  by  an  uncertain  parameter 
vector.  A  single  compensator  that  simultaneously  stabilizes  the  system  is 
sought.  The  conditions  under  which  this  compensator  exists  are  stated  along 
with  the  construction  of  the  compensator.  The  design  is  shown  to  be  suitable 
for  a  class  of  non-minimum  phase  systems  also. 

3.1  TIME  DOMAIN  CONTROL  DESIGN  FOR  ROBUST  STABILIZATION  AND  REGULATION 
There  is  a  considerable  amount  of  literature  on  the  aspect  of  designing 
linear  controllers  for  linear  time  invariant  systems  with  small  parameter 
uncertainty.  However,  for  uncertain  systems  whose  dynamics  are  described  by 
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interval  matrices  (i.e.,  matrices  whose  elements  are  known  to  vary  within  a 
given  bounded  interval) ,  linear  control  design  schemes  that  guarantee 
stability  have  been  relatively  scarce.  Vinkler  and  Wood  [43]  compare  several 
techniques  for  designing  linear  controllers  for  robust  stability  for  a  class 
of  uncertain  linear  systems.  Among  the  methods  considered  are  the  standard 
linear  quadratic  regulator  (LQR)  design,  the  guaranteed  cost  control  (GCC) 
method  of  Chang  and  Peng  [15],  and  the  multistep  guaranteed  cost  control 
(HGCC)  of  Vinkler  and  Vood.  In  these  methods,  the  weighting  on  state  in  a 
quadratic  cost  function  and  the  Riccati  equation  are  modified  in  the  search 
for  an  appropriate  controller.  Also  the  parameter  uncertainty  is  assumed  to 
enter  linearly  and  restrictive  conditions  are  imposed  on  the  bounding  sets. 

In  Chang  and  Peng,  norm  inequalities  on  the  bounding  sets  are  given  for 
stability  but  they  are  conservative  since  they  do  not  take  advantage  of  the 
system  structure.  There  is  no  guarantee  that  a  linear  state  feedback 
controller  exists.  Thorp  and  Barmish  [44]  utilize  the  concept  of  "matching 
conditions"  (MC) ,  which  in  essence  constrain  the  manner  in  which  the 
uncertainty  is  permitted  to  enter  into  the  dynamics  and  show  that  a  linear 
state  feedback  control  that  guarantees  stability  exists  provided  the 
uncertainty  satisfies  matching  conditions.  By  this  method  large  bounding  sets 
produce  large  feedback  gains  but  the  existence  of  a  linear  controller  is 
guaranteed.  But  no  such  guarantee  can  be  given  for  general  "mismatched" 
uncertain  systems.  Hollot  and  Barmish  [45]  and  Schmitendorf  [46]  present 
methods  which  need  the  testing  of  definiteness  of  a  Lyapunov  matrix  obtained 
as  a  function  of  the  uncertain  parameters.  Ackermann  [47],  in  the  multimodel 
theory  approach,  considers  a  discrete  set  of  points  in  the  parameter 
uncertainty  range  to  establish  the  stability.  This  section  addresses  the 
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stabilization  problem  for  a  continuous  range  of  parameters  in  the  uncertain 
parameter  set  (i.e.,  in  the  context  of  interval  matrices).  The  proposed 
approach  attacks  the  stability  of  the  interval  matrix  problem  directly  in  the 
matrix  domain  rather  than  converting  the  interval  matrix  to  interval 
polynomials  and  then  testing  the  Kharitinov  polynomials  (Soh  and  Evans  [48]). 
Towards  this  direction,  bounds  on  the  individual  elements  of  the  perturbation 
matrix  are  developed  in  Yedavalli  [20]  which  are  shown  to  be  less  conservative 
than  the  existing  measures  since  the  method  utilizes  the  structural 
information  of  the  uncertainty,  and  are  applicable  to  systems  with  time 
varying  perturbations.  In  Yedavalli  [49],  a  control  design  method  is 
presented  in  which  the  linear  state  feedback  gain  is  obtained  by  the  standard 
Riccati  equation  and  then  the  robustness  of  the  gain  is  investigated  by  the 
resulting  elemental  perturbation  bound.  In  this  section  we  design  the  gain 
such  that  its  determination  directly  uses  the  uncertainty  structure;  this  is 
done  by  parameter  optimization.  In  other  words,  the  controller  gains  are 
determined  such  that  they  maximize  (in  a  certain  sense)  the  elemental 
perturbation  bound  for  a  given  uncertainty  structure.  The  design  technique  is 
then  extended  to  the  observer  based  feedback  case  as  considered  in  Menga  and 
Dorato  [50]  and  the  relative  tradeoffs  are  discussed. 


1 LfltMfc  ¥  W  jfcl  '4  TH  CO  f  W  iI«U  CO  <  ‘M  M  IH^1] 


In  this  section,  we  first  briefly  review  the  upper  bounds  for  robust 
stability  presented  in  Reference  [20]  for  "structured"  (elemental) 
perturbation.  Structured  perturbations  are  those  for  which  magnitude  bounds 
on  the  individual  matrix  elements  are  known  for  a  given  model  structure. 

Consider  the  following  linear  dynamic  system 
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x(t)  -  A(t)x(t) 


(3.1a) 


-  [A,,  +  E(t)J  x(t)  (3.1b) 

where  x(t)  -*  Rn  is  the  state  vector.  A0  is  the  n  x  n  nominally  stable  matrix 
and  E(t)  is  the  "error"  matrix.  In  the  case  of  structured  perturbation,  the 
elements  of  E(t)  are  such  that 


Max  |  Etj(t)  I  -  c 


+ 1 


.  Max  e. , 

, ,  and  c  -  „  1J 


(3.2) 


In  Reference  [20],  it  is  shown  that  the  system  (3.1)  (with  (3.2))  is 
asymptotically  stable  if 
1 


ciJ  <  a  (P  U  ) 

mu'  m  «/ • 


u 


•u 


-  MU 


•ij 


(3.3a) 


for  all  U#1J  /  0,  i,j  -  l,...,n,  where  P  satisfies  the  Lyapunov  matrix 
equation 

PA0  +  A*P  +  2In  -  0  (3.3b) 


and 


u.ij  *  £ij/£  <Thus  0  <  u.ij  <  X>  (3.3c) 

Here  <?(•)  denotes  the  singular  value  of  (•),  (•),  denotes  the  symmetric  part 
of  the  matrix  (•),  and  (•)„  denotes  the  matrix  formed  with  the  absolute  values 
of  the  entries  of  the  matrix  (•).  Simple  examples  illustrating  this  bound  and 
the  role  of  matrix  U,  in  utilizing  the  structural  information  about  the  error 
matrix  are  given  in  Reference  [20]. 

It  may  be  noted  that  U,  can  be  formed  even  if  one  knows  only  the  ratio 
e1J/£  instead  of  knowing  the  two  terms  separately.  One  suitable  choice  for 
the  ratio  is 

u.U  "  cij/£  -  I Aol j  I /  I A0lj  | max  (3.3d) 

for  all  i,j  for  which  /  0. 
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Remark :  From  (3.2),  it  is  seen  that  ctJ  are  the  maximum  modulus 
deviations  expected  in  the  individual  elements  of  the  nominal  matrix  A0.  If 
we  denote  the  matrix  A  as  the  matrix  formed  with  c1J(  then  clearly  A  is  the 
"majorant"  matrix  of  the  actual  error  matrix  E(t).  It  may  be  noted  that  U#  is 
simply  the  matrix  formed  by  normalizing  the  elements  of  A  (i.e.,  cij)  with 
respect  to  the  maximum  of  c1J  (i.e.,  c). 

i.e.,  A  -  Ut  (absolute  variation)  (3.4) 

Thus  here  are  the  absolute  variations  in  AoiJ.  Alternatively  one  can 
express  A  in  terms  of  percentage  variations  with  respect  to  the  entries  of 
A0lJ.  Then  one  can  write 

A  -  SA^  (relative  or  percentage  variation)  (3.5) 

where  A^j  -  |Aolj  |  for  all  those  i,j  in  which  variation  is  expected  and  6 
are  the  maximum  relative  variations  with  respect  to  the  nominal  value  of 


OiJ 


and 


c  _  Max  c 
‘  -  i,j  5u- 

Clearly,  one  can  then  get  a  bound  on  6  for  robust  stability  as 
1 


S  < 


(3.6) 


m«xl  r  m'i 

where  F  is  the  same  as  in  (3.3b). 


We  now  define,  as  a  measure  of  stability  robustness,  an  index  called 
stability  robustness  index  0^  as  follows: 

Case  A:  LHS  of  (3.3)  or  (3.6)  is  known  (i.e.,  checking  stability  for  given 
perturbation  range).  For  this  case, 

^  6  p  •  £  (or  pt  •  6).  (3.7) 

Thus  ^  >  0  corresponds  to  the  stability  robustness  region. 
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Case  B:  LHS  of  (3.3)  or  (3.6)  is  not  known  (i.e.,  specifying  the  bound).  For 
this  case, 

Psh  b  (or  Mr)  •  (3.8) 

Usins  the  concept  of  state  transformation,  it  has  been  shown  in 
Yedavalli  and  Liang  [41]  that  it  is  possible  to  further  improve  the  bounds  on 
£ .  For  structured  perturbation,  it  is  possible  to  get  higher  bounds  even  with 
the  use  of  a  diagonal  transformation.  This  result  can  be  stated  as  follows. 
Given  the  diagonal  transformation  matrix 
M  -  Diag  [m1,m2,  . . . ,mj 


the  system  of  (3. 


max 

i.j 


m 


U 


•ij 


U 


'•ij 


-  M.U, 


•U 


or 

«  <  K 

where 

K  -  v^„(pnu,). 

and 

PA„  +  A*P  +  21-0 

o  o  n 

and  U#  is  formed  such  that 

°.ij  “  ^  where  "ij  “  l^j  ".U 

and 

£  -  c.j  ,  Aq  -  M-1A0M. 

max 

(‘Tote  that  in  general  0,  f  M'^M.) 

Extension  to  Closed  Loop  Systems 

Consider  the  linear,  time  invariant  system  described  by 
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x  -  Ax  +  Bu 


(3.9) 


x(0)  -  x0 

where  x  is  an  n  x  1  state  vector,  and  the  control  u  is  m  x  1.  The  matrix  pair 
(A,B)  is  assumed  to  be  completely  controllable. 

For  this  case,  the  nominal  closed  loop  system  matrix  is  given  by 
A  -  A  +  BG,  G  -  -R^K/p,.  (3.10a) 

and 

R'1 

KA  +  AtK  -  KB  — BtK  +  Q  -  0,  (3.10b) 

Pc 

and  A  is  asymptotically  stable. 

Here  G  is  the  Riccati  based  control  gain  where  Q,R0  ere  any  given 
weighting  matrices  which  are  symmetric  and  positive  definite  and  pc  is  the 
design  variable. 

The  main  interest  in  determining  G  is  to  keep  the  nominal  closed  loop 
system  stable.  The  reason  the  Riccati  approach  is  used  to  determine  G  is  that 
it  readily  renders  (A  +  BG)  asymptotically  stable  with  the  above  assumption  on 
Q  and  Rq. 

Now  consider  the  perturbed  system  with  linear  time  varying  perturbations 
EA(t)  and  EB(t),  respectively,  in  the  matrices  A  and  B,  i.e., 

x  -  [A  +  EA(t)]x(t)  +  [B  +  EB(t)]u(t)  (3.11) 

Let  AA  and  AB  be  the  perturbation  matrices  formed  by  the  maximum  modulus 
deviations  expected  in  the  individual  elements  of  matrices  A  and  B, 
respectively.  Then  one  can  write 
AA  -  c,Uea  1 

>  (Absolute  variation)  (3.12) 

AB  -  cbUeb  J 
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where  £a  is  the  maximum  of  all  the  elements  in  AA  and  cb  is  the  maximum  of  all 
elements  in  AB.  Then  the  total  perturbation  in  the  linear  closed  loop  system 
matrix  of  (3.10)  with  nominal  control  u  -  Gx  is  given  by 


A  -  AA  +  ABGn  -  t.Ue,  +  £bUebGB  (3.13) 

Assuming  the  ratio  cb/c,  -  £  is  known,  we  can  extend  the  main  result  of  (3.3) 
to  the  linear  state  feedback  control  system  of  (3.9)  and  (3.10)  and  obtain  the 
following  design  observation. 

Design  Observation  1:  The  perturbed  linear  system  is  stable  for  all 
perturbations  bounded  by  ca  and  £b  if 

1 


£*  < 
and  £b  <  cn  where 

P(A  +  BG)  +  (A  +  BG)tP  +  2In  -  0 
Alternatively,  we  can  write 
AA  -  6  A 

(Absolute  variation) 


(3.14a) 


(3.14b) 


(3.15) 


AB  -  *.  Bm 

d  n 

where  Am^  -  |  A±J  |  and  Bm^  -  l®ijl  for  all  i .  j  in  which  variation  is  expected 
and  Amlj  -  0,  Bm^  -  0  for  all  i,j  in  which  there  is  no  variation  expected. 

For  this  situation,  assuming  5b/5,  -  5  is  known,  we  get  the  following  bound  on 
6t  for  robust  stability. 

Design  Observation  2:  The  perturbed  linear  system  is  stable  for  all 
relative  (or  percentage)  perturbations  bounded  by  6a  and  6b  if 


«.  < 


WP»<A>  +  *V.)]. 


and  5b  <  6nt  where  P  satisfies  (3.14b). 


(3.16) 
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Note  that  the  above  expressions  can  be  suitably  modified  if  only  either 


AA  or  AB  is  present. 

Remark :  If  we  suppose  AA  -  0,  AB  -  0  and  expect  some  control  gain 
perturbations  AG,  where  we  can  write 

AG  -  cgUeg  (3.17a) 

then  stability  is  assured  if 


(3.17b) 


a  (P  B  Ue  ) 

In  this  context,  pg  can  be  regarded  as  a  "gain  margin." 

For  a  given  c>1J  and  £blJ,  one  method  of  designing  the  linear  controller 
would  be  to  determine  G  of  (3.10)  by  varying  pc  and  pick  a  control  gain  that 
satisfies  (3.14).  If  et  and  cb  are  not  known  a  priori,  then  one  can  determine 
pc  of  (3.10)  such  that  p  is  maximum.  For  an  aircraft  control  example  which 
utilizes  this  method,  see  Reference  [49]. 


CgntEdL  Design  with  Stability  Robustness  fianaUalnt 


In  the  previous  section,  efforts  were  directed  to  design  a  linear  full 
state  feedback  controller  for  robust  stability.  However,  in  that  treatment, 
the  control  gain  determination  does  not  directly  involve  the  stability 
robustness  criterion  as  a  design  constraint.  Instead,  for  a  predetermined 
linear  control  gain  (obtained  by  many  different  nominal  methods),  the 
perturbation  bound  is  calculated  and  in  the  cases  where  the  parameter 
perturbation  ranges  are  given,  the  stability  robustness  condition  is  checked 
(for  robust  stability).  Even  though  the  bound  p  utilizes  the  structural 
information  of  the  uncertainty,  this  design  procedure  does  not  utilize  the 
structural  information  Ut  in  the  determination  of  the  control  gain  G. 
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In  this  section,  we  attempt  to  solve  the  problem  of  control  design  for 
robust  stabilization  in  a  more  direct  and  general  way  by  formulating  it  as  a 
parameter  optimization  problem.  Instead  of  designing  the  control  gains  by 
nominal  means  and  then  checking  its  stability  robustness  bounds,  we  propose  to 
include  the  stability  robustness  condition  explicitly  in  the  design  procedure 
as  a  performance  measure.  In  this  way  it  is  possible  to  exploit  (in 
principle)  the  uncertainty  structure  U,  in  the  design  procedure. 

Performance  Index  Specification 

An  optimization  problem  to  maximize  the  stability  robustness  bound  p  can 
be  posed  as  follows.  (For  simplicity,  let  us  consider  the  case  AB  -  0.) 
Minimize 
J,  -  a  (P  U  ) 

i.e.,  maximize  (3.18a) 

1 

M  ”  a  (P  U  ) 
subject  to  the  constraints 

P(A  +  BG)  +  (A  +  BG)TP  +  2In  -  0  (3.18b) 

and 

ReAt(A)  -  ReAt(A  +  BG)  <  0  (3.18c) 

where  A(»)  are  the  eigenvalues  of  (•)  and 

u  -  Gx  (3 . 18d) 

We  now  append  the  above  stability  robustness  measure  to  the  standard  quadratic 
performance  index  in  state  x  and  control  u  to  take  the  control  effort 
constraints  into  consideration.  Thus  the  new  performance  index  would  be 
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+  [J  (xtQx  +  utRu)  dtj 
0 


“  £,»«c<p»u.«).  +  V2  trace  KX„;  R  -  ^Rq,  and  Q  >  0 
where  K  satisfies  the  Lyapunov  equation 

K(A  +  BG)  +  (A  +  BG)tK  +  GtRG  +  Q  -  0 


(3.19) 


(3.20a) 


and 

X0  -  x0xj  (3.20b) 
The  optimization  problem  then  is  as  follows:  Find  G  such  that  the  performance 
index 


Ji  -  +  1/2  trace  KXJ 

(3.21a) 

is  minimized  subject  to  the  constraints 

P(A  +  BG)  +  (A  +  BG)tP  +  2In  -  0 

(3.21b) 

K(A  +  BG)  +  (A  +  BG)tK  +  GtRG  +  Q  -  0 

(3.21c) 

and 

ReU^A  +  BG))  <  0 

(3 . 21d 

Modified  Performance  Index 

Note  that  the  above  performance  index  Jx  contains  a  term  involving  the 
maximum  singular  value  as  well  as  a  positive  matrix  Pa.  Even  though  there  are 
algorithms  to  optimize  g  («) ,  optimization  of  an  index  like  the  one  posed 
(P.U.)  is  a  formidable  task  as  it  is  computationally  very  complex.  Hence  we 
intend  to  modify  the  performance  index  such  that  it  becomes  more  tractable. 

Noting  that  the  Frobenius  norm  of  a  matrix  is  always  an  upper  bound  on 
the  spectral  norm  of  the  matrix,  i.e., 

I  l<->  I  I,  >*«<•>.  (3-22) 
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and  that 


*.„<•)  >  W*>.  (3.23) 

we  propose  the  following  upper  bound  Ja  to  be  minimized  instead  of  g,.,(PrU.„), . 

Pr<?P9§J-tL<?n.,3,l: 

Ja  -  1/2  trace  (PWPT  +  PTWPj  >  cr£„(p.u.>.  (3.24) 

for  some  suitable  diagonal  weighting  matrix  W. 

The  diagonal  weighting  matrix  V  is  such  that  Wa  -  0  whenever 

U.ij(j  -  1,2 . n)  -  0  for  a  given  row  i  and  Wlt  -  wt  whenever  UalJ 

( j  -  1 , 2 , . . . , n)  /  0  for  a  given  row  i  and  any  column  j .  Even  though  the 
specification  of  wA  is  crucial  in  establishing  the  upper  bound  property  of  Ja 
as  in  (3.24),  it  turns  out  that  it  is  possible  to  specify  the  wA  >  0  as 
arbitrary  and  transfer  its  implication  in  the  design  to  another  design 
variable,  namely,  pc,  the  weighting  on  the  control  variable. 

Remark:  At  this  stage  of  research,  one  limitation  of  specifying  the  W 
matrix  as  above  is  that  it  reflects  the  uncertainty  structure  (Ua)  only 
partially  in  the  sense  that  is  the  same  (i.e.,  the  same  diagonal  entry) 
irrespective  of  whether  there  are  uncertain  elements  present  in  different  j-th 
locations  or  only  in  one  j-th  location.  However,  for  those  uncertainty 
structures  Ua  which  make  UaUj  diagonal,  we  can  replace 

W  -  aUaUa  (diagonal) (a  is  a  scalar  >  0)  (3.25) 

in  (3.24),  which  then  amounts  to  utilizing  the  structure  of  the  uncertainty 
completely  and  a  acts  as  a  weighting  parameter.  The  forms  of  Ua  which  render 
UaUj  diagonal  are  given  by: 

Case  A:  variations  in  one  row  only 
Case  B:  variations  in  diagonal  elements 
Case  C:  variations  in  antidiagonal  elements 
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Case  D:  no  tvo  varying  elements  are  in  the  same  column.  In  fact,  Cases 
A,  B,  and  C  are  special  cases  of  Case  D. 

Efforts  are  underway  to  prescribe  a  more  versatile  performance  index 
such  that  it  completely  utilizes  the  uncertainty  structure  and  still  maintain 
the  tractability  property. 

We  are  now  in  a  position  to  state  the  problem  of  finding  the  "optimal” 
state  feedback  gain  G  for  robust  stability  and  nominal  regulation  as  follows: 


Minimize 

J  -  {1/2  trace  (PWPT  +  PTWP]  +  1/2  trace  KX<,)  (3.26a) 

subject  to 

P(A  +  BG)  +  (A  +  BG)tP  +  2In  -  0  (3.26b) 

K(A  +  BG)  +  (A  +  BG)tK  +  GjRoG  +  Q  -  0  (3.26c) 

R.At(A  +  BG)  <  0  (3.26d) 

SfllatlanJa:  Parameter  Optimization 


We  approach  the  solution  to  the  above  nonlinear  (quadratic  performance 
index)  programming  problem  by  writing  down  necessary  conditions  and 
investigating  the  solutions  which  satisfy  them.  Using  the  technique  of 
Lagrange  multipliers,  we  transform  the  above  constrained  optimization  problem 
to  an  unconstrained  optimization  problem  by  defining  the  Hamiltonian.  Thus  we 
write  minimize 

(H)  (3.27) 

where  H  is  the  Hamiltonian  given  by 

H  -  trace  (1/2  (PTWP  +  PWPT)  +  1/2  lO^ 

+  L2(PA  +  AtP  +  2In)  +  Lj(KA  +  AtK  +  GTRG  +  Q) )  (3.28) 

and  L1  and  Lj  are  the  Lagrange  multiplier  matrices  corresponding  to  the  two 
matrix  constraints. 
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The  first  order  necessary  conditions  are  given  by 


au 

—  -  (A  +  BG)TP  +  P(A  +  BG)  +  2In  -  0 

o  L« 


3H 

3Lj 


-  (A  +  BG)TK  +  K(A  +  BG)  +  GtRG  +  Q  -  0 


-  (A  +  BG)L*  +  L*(A  +  BG)1  +  PW  +  WP  -  0 
~~  “  (A  +  BG)L*  +  L|(A  +  BG)1  +  Xj  -  0 


~  -  2Bt(PL1  +  KL-.)  +  2RGLj  -  0 


(3.29a) 

(3.29b) 

(3.29c) 

(3.29d) 

(3 . 29e) 


In  arriving  at  these  conditions,  the  matrix  derivative  identities  given  in 
Athans  [SO]  are  used. 

One  can  determine  the  gain  G  by  simultaneously  solving  for  the  above 
equations,  starting  with  an  initial  guess  G0.  Guidelines  for  obtaining 
solutions  to  the  above  type  of  equations  are  given  in  Levine  and  Athans  [51]. 

Remark :  Note  that  when  the  stability  robustness  constraint  is  absent 
(which  is  the  case  by  making  W  ■  0) ,  the  above  problem  formulation  reduces  to 
the  standard  linear  quadratic  regulator  problem  and  the  equations  (3.29)  yield 
the  standard  Riccati  equation  for  the  optimal  control  gain  G. 

However,  with  the  stability  robustness  constraint  present,  the  gain  G  is 
seen  to  be  a  function  of  the  initial  condition  matrix  Xc.  As  pointed  out  In 
Levine  and  Athans  [52],  this  dependence  of  the  controller  on  the  initial 
condition  can  be  removed  by  treating  xo  to  be  a  random  vector  with  zero  mean 
and  uniformly  distributed  over  a  sphere  of  unit  radius,  thereby  considering 
the  worst  case  situation.  Accordingly,  we  can  modify  the  performance  index  as 

J  -  {1/2  trace  [PWPT  +  PTWP]  +  1/2  trace  K)  (3.30) 


AO 


where  it  is  assumed  that  Xe  -  £(x0xB)  -  In  (with  £  being  the  expectation 
operator) . 

Accordingly,  one  of  the  necessary  conditions,  (3.29d),  changes  to 
(A  +  BG)L*  +  L*(A  +  BG)1  +  IB  -  0  (3.31) 

which  then  allows  us  to  express  the  "optimal"  control  given  G  explicitly, 
using  (3.29e)  as 

G  -  -R'1Bt(PLj  +  KL^1  (3.32) 

Extension  to  Observer  Based  Feedback  Controller 

It  may  be  noted  that  the  above  procedure  can  be  readily  extended  to  the 
case  of  observer  based  feedback  control,  with  the  system  description  given  by 
x  -  [A  +  EA(t)]x(t)  +  [B  +  Eg(t) ]u(t)  (3.33a) 

z  -  [M  +  EB(t)]x(t)  (3.33b) 

where  z  is  the  i -vector  of  measurements  and  EA,  Eg,  and  are  the 
"perturbations"  in  the  nominal  matrices  A,  B,  and  M. 

In  this  connection  a  few  remarks  about  the  paper  by  Menga  and  Dorato 
[52]  are  in  order.  The  problem  formulation  of  designing  an  observer  based 
controller  in  this  section  is  similar  in  spirit  to  that  of  Menga  and  Dorato. 
However,  in  that  paper  no  explicit  bounds  on  the  individual  elements  of  the 
perturbation  are  incorporated  as  Is  done  here.  Also,  some  restrictions  are 
placed  on  the  uncertainty  structure  to  fit  it  into  their  proposed  problem 
formulation  (such  as  orthogonality  of  the  uncertainty  matrices,  and  only  a 
single  uncertain  parameter  being  allowed  in  Eg).  The  major  contribution  of 
the  present  work  that  is  significantly  different  from  that  of  Menga  and  Dorato 
is  the  exploitation  of  the  uncertainty  structure  in  obtaining  the  stability 
robustness  bounds.  This  in  turn  results  in  the  consideration  of  two  separate 
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Lyapunov  equations  In  the  problem  formulation  (Equation  (3.26))  as  opposed  to 
only  one  Lyapunov  equation  considered  in  Menga  and  Dorato. 


The  observer  structure  is  the  standard  Luenberger  observer  (Reference 
[53]),  with 

P  -  Fp  +  Hu  +  Dz  (3.34) 

where  p  is  the  estimate  of  x  and  the  matrices  F,  D,  and  H  satisfy  the  observer 
conditions 

SA  -  DM  -  FS  (3.35a) 

H  -  SB  (3.35b) 

for  an  appropriate  transformation  matrix  S.  The  control  is  given  by 

u  -  Gp  (3.36) 

For  brevity,  the  details  of  the  problem  formulation,  which  follows  the 
development  given  in  previous  sections,  are  not  given  here. 

Example  3.3:  Consider  a  simple  second  order  linear  time  invariant 
system  given  by 


0 


a 


1 


-0.5 


x  + 


0 

1 


u 


(3.37) 


where  a  is  the  uncertain  parameter  with  nominal  value  a  -  1.  Notice  that 

0  0 

u..  - 

1  0 


We  select 


It  may  be  noted  that  the  robustness  weighting  matrix  W  incorporates  the 
uncertainty  structure  (that  only  a21  element  is  varying)  in  an  explicit  way. 

Now  with  Q  -  I2,  R  -  pR0  -  p  (i.e.,  R0  -  1  and  p  as  a  design  variable) 
and  a  -  1  for  "robust"  design  and  a  -  0  for  nominal  design,  we  can  get  the 
"optimal"  control  gain  G  by  following  the  proposed  procedure. 

The  comparison  of  "robust  (PO)  state  feedback  control  law"  versus  the 
"nominal  state  feedback  control  law"  is  depicted  in  Figure  1,  where  the 
perturbation  bound  p21  is  plotted  against  the  nominal  control  effort 

■< 

uTudt  1/2 . 

v  0 

As  anticipated  for  a  given  control  effort,  the  robust  control  law  yields  a 
higher  perturbation  bound  nzl  than  the  nominal  control  law,  indicating  the 
usefulness  of  the  proposed  optimization  procedure.  Also,  since  in  this  case 
the  "matching  condition"  is  satisfied,  it  is  seen  that  the  higher  the  control 
effort,  the  greater  is  the  perturbation  bound. 

In  the  next  section,  we  present  a  design  procedure  for  systems 
represented  by  transfer  functions  with  uncertain  coefficients. 

3.2  FREQUENCY  DOMAIN  CONTROL  DESIGN  FOR  ROBUST  STABILIZATION 
(INCLUDING  A  CLASS  OF  NON-MINIMUM  PHASE  SYSTEMS) 

Much  of  the  motivation  for  the  transfer  function  based  frequency  domain 
design  work  comes  from  the  so-called  simultaneous  stabilization  problem 
introduced  in  References  [54]  and  [55].  Given  plants  P0(s),  P1(s) , . . . , Pk(s)  , 


does  there  exist  a  single  compensator  C(s)  that  stabilizes  all  of  them?  In 
their  paper,  Saeks  and  Murray  [54]  develop  geometric  conditions  for 
simultaneous  stabilization  and  state  that  their  solution  is  "mathematical  in 


nature  and  not  intended  for  computational  implementation."  They  indicate  that 
computational  criteria  are  only  known  for  the  two  plant  case  (Reference  [56]). 
The  subsequent  work  of  Vtdyasagar  and  Viswanadham  [55]  is  concerned  with  a 
multi- input  multi-output  (MIMO)  generalization  of  some  of  the  single-input 
single-output  (SISO)  results  of  Reference  [54].  To  this  end,  they  prove  that 
the  problem  of  simultaneously  stabilizing  k+1  plants  is  equivalent  to  the 
problem  of  simultaneously  stabilizing  k  plants  with  the  added  requirement  that 
the  compensator  itself  must  be  stable.  As  far  as  computational  criteria  are 
concerned,  the  results  in  Reference  [55]  imply  a  complete  solution  for  the  two 
plant  case;  i.e.,  upon  reducing  the  two  plant  problem  to  that  of  finding  a 
stable  compensator  for  a  single  plant,  one  can  apply  the  results  of  Youla, 
Bongiorno,  and  Lu  [57],  To  illustrate,  if  P0(s)  and  P1(s)  are  strictly  proper 
SISO  transfer  functions  with  P0(s)  stable,  then  the  results  of  Vidyasagar  and 
Viswanadham  lead  to  the  requirement  that  the  "difference  plant," 

Pz(s)  -  P0(s) ,  be  stabilizable  via  a  stable  compensator.  Hence,  according  to 
Reference  [57],  the  problem  reduces  to  checking  for  satisfaction  of  the  parity 
interlacing  property.  Namely,  we  examine  the  pole  zero  pattern  of 
P1(s)  -  P0(s)  and  require  that  no  zeros  on  the  non-negative  real  axis  lie  to 
the  left  of  an  odd  number  of  real  poles,  multiple  poles  counted  according  to 
their  multiplicity. 

It  is  also  shown  in  Reference  [55]  that  given  two  n  x  m  plants,  one  can 
generically  stabilize  them  simultaneously,  provided  that  either  n  or  m  is 
greater  than  one.  This  result  is  further  generalized  in  Reference  [58],  where 
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it  is  shown  that  general  simultaneous  stabilizability  of  r  n  x  m  plants  is 
guaranteed  if  max(n,m)  >  r. 

In  view  of  the  results  of  Reference  [55] ,  the  issue  of  finding  a 
computationally  feasible  test  for  simultaneous  stabilizability  (for  three  or 
more  plants)  was  raised  again  in  a  paper  by  Emre  [59].  In  his  work,  SISO 
plants  are  considered,  and  the  problem  of  finding  a  computational  test  is 
solved  for  the  special  case  obtained  by  imposing  a  constraint  that  all  k+1 
closed  loop  systems  must  end  up  having  the  same  characteristic  polynomial. 

In  Reference  [60],  Barmish  and  Wei  derive  sufficient  conditions  under 
which  a  family  of  SISO  systems  can  be  simultaneously  stabilized  by  a  proper 
(or  strictly  proper  if  desired)  stable  compensator.  Regularity  conditions  are 
imposed  on  the  plant  family  coefficients,  and  it  is  assumed  that  the  plant 
family  is  minimum  phase  with  one  sign  high  frequency  gain.  A  computation 
procedure  for  constructing  a  robust  compensator  is  also  provided.  These 
results  have  been  further  generalized  for  MIMO  systems  by  Wei  and  Barmish 
[61]. 

In  this  section,  we  intend  to  relax  the  minimum  phase  requirement  and 
show  that  there  exist  certain  classes  of  non-minimum  phase  systems  that  can  be 
simultaneously  stabilized  by  a  single  compensator. 

Preliminary  Definitions  and  Lemmas: 

Definition  3.1 

A  family  of  polynomials  with  different  degrees  is  represented  by  an 
uncertain  polynomial  of  the  form 
d(q) 

f(s,q)  -  £  a1(q)s<,<q)  l,  q  e  Q  (3.38) 

i-0 
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where  Q  is  a  subset  of  Rp,  which  is  an  index  set,  a1(*):  Q  •*  R  are 
coefficient  functions  and  d(»):  Q  -*  (0,1,2,...)  is  the  degree  function. 

An  uncertain  polynomial  f(s,q),  as  in  (3.38),  is  said  to  be  Hurwitz 
invariant  (over  Q)  if  for  each  q  6  Q,  all  the  zeros  of  f(s,q)  lie  in  the 
strict  left  half  plane.  f(s,q)  is  said  to  be  positive  Hurwitz  invariant 
(PHI),  over  Q,  if  f(s,q)  is  Hurwitz  invariant  and  a0(q)  >  0  for  all  q  e  Q. 
Finally,  f(s,q)  is  said  to  be  negative  Hurwitz  invariant  (NHI) ,  over  Q,  if 
-f(s,q)  is  positive  Hurwitz  invariant  (over  Q) . 

Definition  3.2 

An  uncertain  polynomial  f(s,q),  as  in  (3.38)  is  said  to  be  standard  if 
it  satisfies  the  following  conditions: 

(1)  Compact  index  set:  The  index  set  Q  is  a  compact  subset  of  Rp. 

(2)  Bounded  degree:  d(*)  is  bounded  over  Q,  i.e., 

d^  •  max(d(q)  :  q  6  Q)  <  +«. 

(3)  Continuous  coefficients:  Given  any  possible  degree  n  <  d^,  it 
follows  that  a0(»),o1(«) , . . . ,on(»)  are  continuous  over  Q. 

(4)  Closed  degree  sets:  Given  any  possible  degree  n  <  d^,,,,  it  follows 
that  the  degree  set 

Q„  -  (q  e  Q:  d(q)  -  n) 

is  a  closed  subset  of  Q.  In  view  of  the  fact  that  Q  is  assumed  to  be  compact, 
closedness  of  Qn  implies  that  Qn  is  compact;  i.e.,  a  closed  subset  of  a 
compact  set  is  compact. 

Remark:  The  condition  of  closed  degree  sets  has  the  effect  of 
partitioning  a  family  of  polynomials  into  closed  subsets  each  having  its  own 
distinct  degree.  A  simple  example  which  violates  the  condition  is 

f(s,q)  -  qs2  +  s  +  1;  q  e  [0,1]. 
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With  Q  -  [0,1],  we  obtain  Qj  -  [0],  but  Q2  -  (0,1],  which  is  not  closed. 


Definition  3.3 

An  uncertain  proper  rational  function  is  defined  to  be  a  ratio  of  the 

form 

P(s,q)  ■  f1(s ,q)/f2(s ,q) ,  q€Q  (3.39) 

where  f1(s,q)  and  f2(s,q)  are  uncertain  polynomials  as  in  (3.38)  and  for  each 
q  e  Q,  deg[Fj(s,q)]  <  deg[ f2(s ,q) ] .  An  uncertain  proper  rational  function,  as 
in  (3.39),  is  said  to  be  standard  if  both  f^s.q)  and  f2(s,q)  are  standard 
uncertain  polynomials. 

Lemma  3.1  (see  Reference  [62]  for  proof) 

Assume  f(s,q),  as  in  (3.38),  is  an  n-th  order  standard  PHI  polynomial 
and  f ' (s , q)  is  an  n'-th  order  standard  uncertain  polynomial  satisfying  the 
following  conditions: 

( 1 )  n'  <  n  +  1  and 

(2)  if  n'  -  n  +  1,  then  <*g(q)  >  0  for  all  q  6  Q. 

Then,  there  exists  an  c*  >  0  such  that  for  any  c  e  (0,c*), 

fc(s,q)  -  f (s , q)  +  cf'(s,q) 
is  also  PHI  over  Q. 

The  above  lemmas  can  be  easily  generalized  for  uncertain  polynomials 
with  different  degrees. 

Lemma  3 . 2  (See  [62]  for  proof.) 

Assume  f(s,q),  as  in  (3.38),  is  a  n-th  order  standard  PHI  polynomial, 
and  f'(s,q)  is  a  n'-th  order  standard  uncertain  polynomial  satisfying  the 
following  conditions: 

(I )  n '  <  n  +  1 ,  and 

(2)  an’ (s > q)  >  0  for  a11  q  €  q. 
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Then,  there  exists  an  c*  >  0  such  that  for  any  c  €  (0 ,c*), 
ft(s,q)  -  sf(s,q)  +  cf'(s.q) 
is  also  PHI  over  Q. 

I-gflima  3.3  (See  Appendix  D  for  proof.) 

Assume  f(s,q),  as  in  (3.38),  is  a  standard  PHI  polynomial,  and  f'(s,q) 
is  a  standard  uncertain  polynomial  satisfying  the  following  conditions:  for 
each  q  e  Q, 

(1)  d'(q)  <  d(q)  +1,  and 

(2)  if  d' (q)  -  d(q)  +  1  for  some  q  6  Q,  then  a{,(q)  >  0. 

Then,  there  exists  an  c*  >  0  such  that  for  any  c  e  (0,c*), 

f4(s,q)  -  f (s , q)  +  cf'(s.q) 
is  also  PHI  over  Q. 

Lemma  3.4  (See  Appendix  D  for  proof.) 

Assume  f(s,q),  as  in  (3.38),  is  a  standard  PHI  polynomial  and  f'(s,q)  is 
a  standard  uncertain  polynomial  satisfying  the  following  conditions: 

( 1 )  d' (q)  <  d(q)  +  1  and 

(2)  e»j,(q)(q)  >  0  for  all  q  e  Q. 

Then,  there  exists  an  c*  >  0  such  that  for  any  e  €  (0,£*), 
ft(s,q)  -  sf (s,q)  +  cf'(s,q) 
is  also  PHI  over  Q. 

Robust  Stabilizabilltv  of  SISO  Uncertain  Systems 

Consider  a  SISO  uncertain  system  represented  by  a  family  of  transfer 
functions  (P(s,q)}:  qeQ: 
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P(s,q) 


(3.40) 


(1)  Standardness:  P(s,q)  is  a  standard  uncertain  rational  function. 

(2)  One  sign  high  frequency  gain:  The  sign  of  the  ratio  of  the  leading 
coefficients  of  D(s,q)  and  N(s,q)  is  sign  invariant  over  Q.  In  the  sequel,  we 
keep  the  sign  of  the  leading  coefficient  of  D(s,q)  positive. 
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(3)  Minimum  phase:  For  each  fixed  q  €  Q,  P(s,q)  is  of  minimum  phase; 
i.e.,  all  zeros  of  P(s,q)  lie  in  the  strict  left  half  plane. 

Theorem  3 ■ 1  (see  Reference  [60]  for  proof) 

Under  the  Assumptions  3.1,  a  SISO  uncertain  system  P(s,q)  is 
simultaneously  stabilizable.  Furthermore,  a  robust  stabilizer  C(s)  can  be 
constructed  to  be  proper  (or  strictly  proper,  if  desired)  and  stable. 
Assumption  3.2 

Assume  that  the  transfer  function  P(s,q)  satisfies  the  following 
conditions : 

(1)  Standardness:  P(s,q)  is  a  standard  uncertain  rational  function. 

(2)  One  sign  high  frequency  gain):  The  sign  of  the  ratio  of  the  leading 
coefficients  of  D(s,q)  and  N(s,q)  is  sign  Invariant  over  Q. 

(3)  The  numerator  polynomial  N(s,q)  -  sNx(s,q)  and  Nx(8,q)  is  Hurwitz 
invariant  over  Q;  i.e.,  all  zeros  of  N^s.q)  lie  in  the  strict  left  half 
plane. 

(A)  “o(q)  “i-(q)(q)  >  0  for  all  q  6  Q. 

Theorem  3 . 2  (See  Appendix  E  for  proof.) 

Under  the  Assumptions  3.2,  a  SISO  uncertain  system  P(s,q)  is 
simultaneously  stabilizable.  Furthermore,  a  robust  stabilizer  C(s)  can  be 
constructed  to  be  proper  (or  strictly  proper,  if  desired)  and  stable. 
Assumption  3.3 

Assume  that  the  transfer  function  P(s,q)  satisfies  the  following 
conditions: 

(1)  Standardness:  P(s,q)  is  a  standard  uncertain  rational  function. 

(2)  One  sign  high  frequency  gain:  The  sign  of  the  ratio  of  the  leading 
coefficients  of  D(s,q)  and  N(s,q)  is  sign  invariant  over  Q. 
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(3)  The  denominator  polynomial  D(s,q)  -  sDx(s,q)  and  Dx(s,q)  is  Hurwitz 
invariant  over  Q,  i.e.,  all  zeros  of  Dx(s,q)  lie  in  the  strict  left  half 
plane . 

(4)  Q0(q)a1B(q)(q)  >  0  for  all  q  €  Q. 

Theorem  3.3  (See  Appendix  F  for  proof.) 

Under  the  Assumptions  3.3,  a  SISO  uncertain  system  P(s,q)  is 
simultaneously  stabilizable.  Furthermore,  a  robust  stabilizer  C(s)  can  be 
constructed  to  be  proper  (or  strictly  proper,  if  desired)  and  stable. 
Assumption  3.4 

Assume  that  the  transfer  function  P(s,q)  satisfies  the  condition  that 
P(s ,q)  -  P0(s)  +  Pl(s,q),  where  P„(s)  is  stable  and  P1(s,q)  satisfies  either 
Assumptions  3.1,  3.2,  or  3.3,  respectively . 

Theorem  3.4  (See  Appendix  G  for  proof.) 

Under  the  Assumptions  3.4,  a  SISO  uncertain  system  P(s,q)  is 
simultaneously  stabilizable  by  a  proper  (or  strictly  proper,  if  desired) 
compensator . 

The  Stabilizing  Compensator  Synthesis  and  Illustrative  Example 

In  this  section,  we  extract  the  crucial  Ingredients  for  compensator 
construction  from  the  proofs  of  Theorem  3.2,  3.3,  and  3.4. 

Procedure  of  Synthesis  3.1 

Assume  that  the  system  satisfies  Assumption  3.1.  The  compensator 
computation  procedure  is  given  in  Reference  [60]. 

Procedure  of  Synthesis  3.2 

Assume  that  the  system  satisfies  Assumption  3.2.  Design  compensator 


Cx(s)  -  Nc(s)/De  0(s)  by  following  Procedure  3.1  such  that 
A0(s,q)  ■  Nff(s)Nj(s , q)  +  De  0(s)D(s,q) 
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is  PHI  over  Q.  Given  A1(s,q),  select  et  >  0  such  that 
Ain(s,  q)  -  sA^s.q)  +  e1Dc.0(s)D(s,q) 
is  PHI  and 

Dc.ih(s)  -  *osi  +  cisW  +  ...  +  Ci 

is  Hurwitz.  Assumption  3.2  guarantees  the  success  of  this  procedure  for 

i  -  0,1,2 . A-l.  The  compensator  is  given  by 

C(s)  -  Nc(s)/Dc  A(s)De  0(s). 

Procedure  of  Synthesis  3.3 

Assume  that  the  system  satisfies  Assumption  3.3.  Set 
A0(s,q)  ■  D1(s,q) 

which  is  PHI  over  Q.  Given  A1(s,q),  select  51  >  0  such  that 
di+1(s,q)  ■  sAt(s , q)  +  ^(s.q) 

is  PHI.  Assumption  3.3  guarantees  the  success  of  this  procedure  for 

i  -  0,1,2 . A-l.  The  numerator  of  the  compensator  is  given  by 

Nc(s)  -  50s'1  +  5jS'2  +  ...  +  8X_ j. 

Set 

AJ(s)  -  Ne(s)N(s,q)  +  sAD(s,q) 

For  given  Ai(s,q),  select  et  >  0  such  that 

+  e^tys.q) 

is  PHI  and 

Dc.i<s>  ■  q  +  +  . . .  +  «4sl 

is  Hurwitz.  Assumption  3.3  guarantees  the  success  of  this  procedure  for 
i  -  1,2,. ..,1  >  A-l.  The  compensator  is  given  by 
C(s)  -  Ne(s)/De ^(s). 
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Procedure  of  Synthesis  3.4 

Assume  that  the  system  satisfies  Assumption  3.4.  Design  a  compensator 
Cx(s)  -  Nc  1(s)/Dc  1(s)  such  that 

A0(s,q)  -  Ne^  x(s)Nx(s ,  q)  +  De  <1(s)D1(s  ,  q) 
is  PHI  over  Q.  Then  the  desired  compensator  C(s)  is  given  by 


C(s)  - 


C,(s) 


1  -  C1(s)P0(s) 

In  the  computation  procedures,  it  is  often  required  to  determine  Hurwitz 
invariance  of  uncertain  polynomials.  Numerous  methods  can  be  found  in 
References  [2]  through  [5],  as  well  as  the  results  of  subsection  2.2. 


A  continuum  of  plants  is  described  by 

P(s,q)  -  P0(s)  +  P1(s,q)  -  +  2  ^  ,  qx  6  [1,2],  q2  €  [-2,2]. 

S  T  1  S  *  498  '  A, 

Since  Assumption  3.4  is  observed  to  be  satisfied  by  inspection,  this  family  of 
plants  is  simultaneously  stabilizable .  Ve  now  construct  a  stabilizing 
compensator.  Design  C0(s)  -  N„(s)/D0(s)  such  that 

A0(s,q)  -  N0(s)(s  +  qx)  +  (D„(s)(s2  -  q2s  +  1)) 
is  Hurwitz  invariant.  Following  Procedure  3.1,  we  take  N0(s)  -  1  and 
D0(s)  -  c0;  then 


A0(s,q)  -  (s  +  qx)  +  £0(s2  -  q2s  +  1)  -  £0s2  +  (1  -  £0q2)s  +  £„  +  qx 
Obviously,  c0  <  0.5  will  make  A0(s,q)  Hurwitz  invariant,  so  we  choose 
c q  -  0.4.  Now  consider 

A1(s,q)  -  sA0(s,q)  +  £x(s2  -  q2s  +  1) 
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It  is  easy  to  check  (using  the  Kharitonov  test  [2],  for  instance)  that  when 


£x  -  0.4,  A1(s,q)  will  be  Hurvitz  invariant.  Hence,  the  robust  compensator 
for  P1(s,q)  is 


Cx(s) 


1 _ 2.5 

0.4s  +0.4  s+1 


Consequently,  the  robust  compensator  C(s)  for  P(s,q)  is  given  by 


C(s) 


Ct(s) 

1  -  Cx(s)P0(s) 


2. 5(s  +  1) 
s2  +  2s  -  1.5 
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IV.  ROBUST  CONTROL  OF  LINEAR  SYSTEMS  WITH  COMBINED  STRUCTURED  (PARAMETRIC) 
AND  UNSTRUCTURED  (UNMODELED  DYNAMICS)  UNCERTAINTY 

Most  of  the  published  literature  on  the  robust  stabilization  (and 
performance)  problem  treats  the  two  cases  of  structured  uncertainty 
(parametric  variations)  and  unstructured  uncertainty  (unmodeled  dynamics) 
separately.  The  aspect  of  unstructured  uncertainty  is  conveniently  handled  in 
the  frequency  domain  framework,  whereas  the  structured  uncertainty  problem  is 
treated,  as  discussed  in  the  previous  chapters,  both  in  frequency  domain  and 
time  domain  frameworks .  While  research  in  each  of  these  categories  is  still 
growing  at  an  impressive  pace,  with  a  variety  of  problems  being  addressed  and 
resolved,  the  motivation  for  carrying  out  research  that  treats  the  combined 
case  in  a  unifying  framework  comes  from  a  variety  of  reasons.  Many 
engineering  applications  naturally  involve  the  presence  of  both  types  of 
uncertainties;  the  application  of  large  space  structure  control  is  a  case  in 
point.  Also,  the  mathematical  frameworks  considered  for  these  two  cases  of 
uncertainty  cannot  always  be  transformed  from  one  to  the  other.  Hence,  in 
this  final  leg  of  the  research  work  done  under  this  contract,  we  consider  the 
combined  uncertainty  problem.  In  Section  4.1,  we  concentrate  on  the  Large 
Space  Structure  models  in  the  state  space  framework  and  model  the  unmodeled 
dynamics  error  as  an  additive  perturbation,  and  using  the  theory  developed  by 
the  author,  obtain  bounds  on  the  control  gains  to  keep  the  closed  loop  system 
stable  under  both  control  and  observation  spillover.  Then  the  aspect  of 
parameter  variations  is  addressed  and,  using  the  same  tool  as  before,  a  robust 
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controller  is  designed  for  stabilization  under  these  parameter  variations. 

Then  in  section  4.2  a  frequency  domain  based  design  method  is  given  to 
construct  a  single  compensator  for  robust  stabilization  in  the  presence  of 
combined  structured  and  unstructured  uncertainties.  Apart  from  the  difference 
in  perspectives,  another  main  difference  between  these  two  methods  is  that  in 
the  time  domain  design,  the  two  types  of  uncertainties  are  accommodated  in 
sequence,  whereas  in  the  frequency  domain  method,  simultaneous  presence  is 
addressed. 


4.1  TIME  DOMAIN  (STATE  SPACE)  DESIGN  METHOD  FOR  LSS  MODELS 

One  fundamental  design  problem  in  the  control  of  LSS  is  the  control  of 
large  dimensional  system  with  a  controller  of  much  smaller  dimension.  In  this 
framework,  the  modeling  error  is  "mode  truncation."  One  prominent  effect  of 
mode  truncation  in  the  control  design  process  is  the  instability  caused  due  to 
the  interaction  with  residual  modes,  labeled  "spillover"  [63].  Furthermore, 
the  parameters  of  the  control  design  mod?l  Itself  (such  as  modal  frequencies, 
mode  shapes,  and  damping  ratios)  are  all  known  imprecisely.  For  instance,  the 
errors  associated  with  the  finite  element  models  of  LSS  Increase  typically 
with  Increase  in  the  mode  number  (l.e.,  tiui  muuei  frequency).  Even  though  in 
most  cases  the  effect  of  this  "parameter  variation”  error  is  performance 
degradation,  it  is  sometimes  possible  to  encounter  instability  also.  In  the 
light  of  these  observations  it  is  evident  that  "robustness"  is  an  extremely 
desirable  feature  of  any  active  control  design  proposed  for  LSS  control. 

The  bulk  of  the  published  literature  on  the  "robust  stability"  of  LSS 
models  addresses  the  modeling  error  "mode  truncation"  (i.e.,  the  use  of 
reduced  order  models/controllers)  in  the  design  [64] -[68].  The  problem  of 
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"parameter  errors"  in  the  control  design  of  LSS  has  not  received  as  much 
attention  as  the  spillover  problem.  In  the  former  category,  recent  results 
using  the  frequency  domain  viewpoint  treat  the  mode  truncation  error  as  an 
unstructured  (high  frequency)  uncertainty  and  propose  solution  methods  such  as 
LQG/LTR  [69] -[70].  Applications  of  LQG/LTR  method  in  LSS  control  problems  are 
discussed  in  References  [71]  and  [72].  In  the  time  domain  viewpoint,  Ikeda 
and  Siljak  [73]  give  bounds  on  the  spillover  terms  using  vector  Lyapunov 
functions.  Time  domain  control  design  schemes  which  treat  both  mode 
truncation  and  parameter  variations  for  LSS  control  have  been  relatively 
scarce  [74]-[75],  and  even  among  these,  Reference  [74]  treats  only  the  small 
parameter  variation  case  while  Reference  [75]  assumes  stability  and  addresses 

*• 

the  performance  aspect.  The  proposed  method  is  a  contribution  in  this 
direction.  In  this  section  both  the  "mode  truncation  error"  and  the 
"parameter  error"  are  treated  as  an  additive  perturbation,  and  the  stability 
robustness  analysis  for  structured  uncertainty  developed  by  the  author  in 
References  [77]  and  [20]  is  extended  to  treat  the  above  problem  to  obtain 
bounds  on  the  "spillover"  terms  to  maintain  stability.  In  this  context,  the 
methodology  of  this  paper  is  similar  in  spirit  to  the  paper  by  Ikeda  and 
Siljak,  but  with  considerable  difference  in  detail.  One  important  feature  of 
the  proposed  method  is  that  the  bounds  on  control/estimator  gains  can  be 
determined  apriori,  simply  using  the  open  loop  modal  data;  also,  the  tradeoff 
between  these  bounds  and  the  number  of  modes  one  wishes  to  control  is  clearly 
brought  out. 

Description  of  the  Proposed  Approach 

As  noted  in  the  previous  section,  the  two  major  design  issues  of  concern 
in  LSS  control  are  those  of  the  avoidance  of  spillover  [65],  and  parameter 
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uncertainty.  Spillover  is  brought  about  by  a  "two-model"  design  philosophy  in 
which  a  higher  order  "evaluation  model"  (which  serves  the  role  of  the  physical 
system  in  system  evaluation)  is  driven  by  a  lower  order  controller  of 
sufficiently  large  control  effort.  Parameter  uncertainty  arises  because  of 
the  inherent  inaccuracies  in  the  development  of  mathematical  models  (e.g. 
finite  element  modeling).  In  the  proposed  design  scheme,  both  of  these  issues 
are  addressed  in  a  systematic  way.  First  the  aspect  of  mode  truncation  error 
(spillover)  is  addressed.  While  there  are  many  methods  available  for 
obtaining  the  reduced  order  controller  [75] -[76],  we  employ  the  route  of 
reducing  the  open  loop  evaluation  model  (by  some  known  criterion)  and  then 
building  a  full  order  linear  quadratic  gaussian  (LQG)  controller  for  the 
reduced  order  model.  We  then  apply  the  results  of  the  Elemental  Perturbation 
Bound  Analysis  (EPBA)  discussed  in  section  3.1  to  the  closed  loop  system  model 
for  obtaining  upper  bounds  on  the  control/estimator  gains  as  a  function  of  the 
number  of  modes  retained  in  the  reduced  order  control  design  model  to 
guarantee  avoidance  of  spillover  [78].  Then,  armed  with  these  bounds  on  the 
gains,  at  the  control  design  stage  (on  the  reduced  order  model),  once  again 
the  results  of  EPBA  are  used  to  obtain  upper  bounds  on  the  tolerable 
variations  in  the  control  design  model  parameters  to  maintain  stability  [79]. 
Thus  the  proposed  "robust  controller"  is  that  controller  which  possesses  as 
high  a  tolerable  perturbation  in  the  system  parameters  as  possible  along  with 
the  guarantee  of  avoidance  of  spillover  instability.  The  important  features 
of  the  proposed  controller  are  that(l)  the  control/estimator  gain  bounds  are 
calculated,  apriori,  using  open  loop  modal  data;  (2)  that  it  Is  a  practically 
implementable  estimator  based  control  law;  and  most  importantly,  (3)  that  it 
brings  about  the  tradeoff  between  the  nominal  performance  attainable  and  the 
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number  of  modes  included  in  the  control  design  model  (i.e.,  the  dimension  of 
the  reduced  order  model) . 

LSS  Models  and  Mode  Truncation  as  Perturbation: 

Consider  the  standard  state  space  description  of  an  LSS  evaluation  model 
with  N  elastic  modes: 


x  -  Ax  +  Bu  +  Dw 
y  -  Cx 
z  -  Mx  +  v 


x(0)  -  x0;  x  -*  R“"  ,  u  -»  R1 
y  -»  Rk 


z  -*  R 


where 


(4.1a) 


x1  -  [x*,x*,  .  .  .  ,xj] ;  xA  - 


*1 
1 1 


(4.1b> 


A  -  Block  diag.  [ .  .  .A...  .  .  ]  ,  A,.  - 


0  1 
'“>i  -2fw1 


(4.1c) 


B1  -  [B^B* . Bj];  B1  - 


(4. Id) 


DT  -  (Dj.D* . Dj] ;  D4  - 


dl 


(4. le) 


C  -  [C1>C2 . C,]  and  M  -  [Mj.Mj . M*] 


(4. If) 


w(t) 

V(t) 


[wI(r)vI(r)] 


tf 

0 


0 

*.V  o 


?(t  -  r) 


(4 . lg) 
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Here  y  is  the  vector  of  variables  we  wish  to  control,  z  is  the 
measurement  vector,  w  and  u  are  zero  mean  white  noise  processes  with  constant 


covariances  W  and  p#V0  respectively,  where  p,  >  0  is  taken  as  a  design 
variable  and  £  is  the  expectation  operator  and  S  is  the  Dirac  delta  function. 

Assuming  vibration  suppression  of  the  flexible  structure  to  be  the 
control  objective,  a  quadratic  performance  index  may  be  written  as 


t  N 

J  -  lim  i  J  wjflf  +  \\  +  PeuTu]  dt 


(4.2) 


which  can  be  written  in  the  form 


J  -  lim  -p  |  f  (yTy  +  peuTu)  dt  -  lim  7  f  [  (yTy  +  pcuTu)  dt 

C  J  ^  J 


-  Jy  +  PCJ 


where  the  matrix  C  of  (4.1)  is  given  by 
C  -  Block  diag.  [...Ci...] 


(4.3) 


(4.4a) 


wi  0 


0  1 


(4.4b) 


The  scalar  pc  >  0  is  also  taken  as  a  design  variable. 

Assuming  the  modes  ija  are  ordered  in  increasing  order  of  frequency,  let 
us  retain  the  first  Nc  modes  for  control  design  purposes.  Accordingly,  the 
reduced  order  control  design  model  of  dimension  ng(-  2Nr)  <  n  is  given  by 
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xR  -  ArXr  +  BrU  +  Drw  Xr  -  R  R  (4.5) 

yR  “  CrXr  zr  “  MrXr  +  v 

where  the  above  control  design  model  is  obtained  by  direct  truncation  of  the 
full  order  model  given  by 


x  - 


Xr 

Ar  o 

Xr 

h 

°R 

- 

+ 

u  + 

Xt 

1 

o 

4* 

Xt 

Bt  _ 

DT  _ 

(4.6) 


y  -  [Cr  ct] 


*R 

>4 


z  -  [Mr  M*] 


*R 


Let  the  full  order  optimal  control  for  the  reduced  order  model  be  given 


by 


u  -  gr  Xr 

xR  “  Ar^r  +  BrU  +  Fr(z  -  MrJCr)  (4.7) 

“  ArXr  +  Frz  where  \  “  A*  +  Brgr  -  FA 
where  GR  and  FR  are  the  standard  "controller"  and  "estimator"  gain  matrices 
obtained  by  minimizing  the  performance  index 

t 


-  linA  ?  [  [  I  wJijJ  +  ]  +  PcuTu"j  dt 

t-~>  {  L  M-l  J  J 


(4.8a) 


-  lim  ~  i  J  (x^cJCrXr  +  pcuTu)  dt  -  JyR  +  pe Ju 
t-*«  i 


(4.8b) 


and  are  given  by 
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(4.9a) 


*A  +  A& 


t,. 


*  <?C,  -  o 


and 


(4.9b) 


Pc 


Txt-1 


Fr  -  —  Wo 


r-1 


KA  +  -  K.M5  -f-  MjjK,  +  D*VD*  -  0 

The  closed  loop  system  matrix  for  the  control  design  model  given  by 


Ar  “ 


(4.9c) 

(4 . 9d) 


Ar 

®R^R 

FrMr 

Ar 

(4.10) 


is  asymptotically  stable  under  the  usual  assumptions  of  controllability  and 
observability. 

The  closed  loop  system  for  the  evaluation  model  is  obtained  by  forcing 
the  evaluation  model  with  the  controller  of  the  control  design  model.  Thus, 
we  have 


*R 

Ar  BrGr  0 

• 

XR 

mm 

FrHr  FrUj 

• 

*T 

0  BtGr  Aj 

y 

[  =.  ®  ct 

u 

o 

cf 

o 

Dr  0 

W 

+ 

o  fr 

*T 

dt  0 

V 

(4.11a) 

[  **  **  *T  l1 


i.e.  , 
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(4.11b) 


"  W  +  ^bw>  Yb  “  ^b^b 
where  the  stability  of  the  matrix  \  is  dictated  by  the  spillover  terms  BTGR 

and  FrMj. 

We  now  write 


\ 

0 

0 

0 

Vs. 

0 

\  -  A  +  E,  - 

0 

A* 

0 

+ 

BrGr  - 

FrMt 

0 

0 

At 

0 

btgr 

0 

Note  that  the  Ej,  matrix  contains  the  spillover  terms  FgHj  and  BTGR  and  that  the 
A  matrix  is  asymptotically  stable  with  block  diagonal  structure.  In  other 
words,  the  spillover  terms  are  treated  as  an  additive  perturbation.  Let  us 
denote  the  matrix  Ej,  as  the  "spillover  matrix." 

It  is  to  be  noted  that  the  matrix  Eb  is  not  known  prior  to  the  control 
design  and  that  A  is  essentially  the  open  loop  plant  matrix. 

Let 

f  -  Max(  |FR1J|),  g  -  Max(  Ig^j  I)  (4.13) 

b  -  Max(  |BU  |  )  ,  m-Max(|Mlj|)  (4-14) 

and 

h  -  Max(f , g) 

Then  the  maximum  element  (in  the  worst  case  situation)  in  the  matrix  Eb  will 
be 

c  -  h(mb  +  im)  -he  (4.15) 

where  the  constant  c  can  be  computed,  apriori,  from  the  given  model.  We  thus 
write  the  "majorant"  of  Ej,  as 

A,,  -  cU  (4.16) 
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where 


U 


0  u#1  0 

U.2  U.3  U.* 


0 


u 


•3 


(4.17) 


with  UBllJ  -  1  whenever  there  is  a  nonzero  entry  in  the  B^  matrix  and 
UBliJ  -  0  for  all  zero  entries  and  so  on  for  the  other  matrices  U>2,  .  .  .  ,U#5. 

Now  using  the  robust  stability  result  of  subsection  3.1  (specifically, 
equation  (3.3)),  we  conclude  that  the  closed  loop  system  matrix  Aj,  is  stable 
if 


£  -  h  c  < 


where  P  satisfies 

PA  +  AtP  +  21  -  0 


% 


(4. 18a^ 


(4.18b) 


Note  that  the  above  equation  can  be  easily  solved  exploiting  the  structure  of 
the  A  matrix.  Accordingly, 

P  -  Block  diag.  [...P41...]  (4.18c) 


l_f_4£2 
2f«i  +  2f 


1/«J 


1  +  u>\ 


(4. 18d) 


Note  that  in  this  case  Plla  -  Pu  (i.e.,  all  the  elements  of  Pu  are  positive). 

Denoting  ^  as  the  bound  obtained  as  in  (4.18a)  when  contains  a 
modes  (a  -  1,2, . . . ,N-1) ,  we  can  get  the  bound  on  the  control/estimator  gain  as 
a  function  of  the  number  of  controlled  modes,  as 
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(4.19) 


k  <  -  >v 

Thus  any  controller/estimator  gains  (G„  and  F„  where  G,  is  n  x  2a  and  FR  is 
2a  x  i)  such  that 

KJ  <*W  |F„.l<h_  <4-20) 

would  keep  the  closed  loop  system  matrix  stable  (i.e.,  avoid  spillover 
instability) . 

Remark :  The  most  important  point  to  note  is  that  the  bound  h^  is 
obtained  apriori  before  control  design  by  simply  utilizing  the  open  loop  modal 
data  such  as  the  modal  frequency,  mode  shape  slopes  at  actuator/sensor 
locations  (B  and  M  matrices),  open  loop  modal  damping  ratio  6,  the  number  of 
control  inputs,  the  number  of  measurements,  and  the  number  of  modes  included 
in  the  control  design  model. 

Since  with  gains  satisfying  (4.20)  (obtained  by  design  variables  pc  and 
p#)  the  closed  loop  system  matrix  Aj,  is  stable,  one  can  proceed  to  compute  the 
nominal  output  regulation  cost  Jy  and  the  corresponding  control  effort  Ju  and 

obtain  a  plot  of  Jy  vs.  Ju  for  each  a  -  1,2 . N-l. 

Stability  of  the  Control  Design  Model  under  Parameter  Uncertainty 

Once  the  stability  under  mode  truncation  is  taken  care  of,  the  next 
issue  is  the  stability  of  the  control  design  model  itself.  From  the  control 
design  equations  of  (4 . 8) - (4 . 20) ,  it  is  seen  that  under  the  nominal  situation 
the  closed  loop  reduced  order  control  design  model  is  asymptotically  stable. 
However,  in  an  LSS  control  problem  the  parameters  of  the  plant  matrix  AR, 
namely,  the  modal  frequencies  and  dampings  as  well  as  the  mode  shape  slopes  at 
actuator  (Br)  and  sensor  (Mr)  locations  are  known  to  be  uncertain.  It  is  also 
known  that  the  uncertainty  in  these  parameters  tends  to  increase  with  an 
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increase  in  the  mode  number.  Thus  with  variations  AA*,  ABg,  AM^  etc.,  the 
control/estimator  set  obtained  in  the  previous  section  cannot  guarantee 
stability  of  the  reduced  order  closed  loop  system. 

Clearly  the  appropriate  step  then  is  to  pick  that  set  of 
control/estimator  gains  which  not  only  satisfy  condition  (4.20),  but  also 
maximize  the  allowable  range  of  variations  in  AA*,  ABg,  and  AMr  for  robust 
stability.  Towards  this  direction,  the  stability  robustness  index  ^SR  of 
(3.7)  or  (3.8)  can  be  used  as  a  robustness  measure  to  be  maximized  with 
parameter  uncertainty  treated  as  an  additive  perturbation.  In  doing  so,  one 
can  even  utilize  the  uncertainty  structure,  namely  the  fact  that  the 
uncertainty  in  modal  parameters  increases  with  increase  in  mode  number.  One 
way  of  characterizing  this  structure  is  as  follows  (here  we  employ  the 
relative  variation  format  of  (3.5)): 
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(And  similarly  for  AM^) 

where  indicate  the  nominal  entries  corresponding  to  the  i-th  mode.  Ue 
assume  -  $b,  which  are  not  known. 

With  the  above  uncertainty  structure,  again  using  the  result  for  robust 
stability  of  (3.3),  one  can  get  a  bound  on  6 a  as 

«.  <  (4.22) 

where  n%  is  obtained  as  usual  by  solving  a  Lyapunov  equation  of  the  type  (4.5) 
using  the  nominally  asymptotically  stable  closed  loop  system  matrix  AR.  Note 
that  ut  is  a  function  of  the  control/estimator  gains.  We  can  thus  define 
as  the  robustness  measure  for  parameter  uncertainty,  i.e., 

4a  *  (4.231 

and  plot  0^  vs.  the  control  effort  Ju. 

Finally  the  proposed  robust  controller,  for  a  given  a  -  1,2,...,N-1  is 
that  controller  which  satisfies  (1)  condition  (4.29),  (2)  Jy  <  (maximum 
tolerable  output  regulation  cost),  and  (3)  as  high  0^  as  possible. 

Of  course  if  there  is  no  stringent  requirement  on  Jy,  then  one  can 
select  a  gain  set  that  achieves  maximum  0^. 

4.2  FREQUENCY  DOMAIN  (TRANSFER  FUNCTION)  DESIGN  METHOD 

In  this  section,  we  consider  the  simultaneous  presence  of  plant 
uncertainties  characterized  by  low  order  parameter  (structured)  variation  and 
high  order  (unstructured)  uncertainty  for  systems  described  by  transfer 
functions .  When  the  system  uncertainty  is  mainly  caused  by  high  order 
unstructured  uncertainty,  the  most  popular  design  model  is  characterized  by  a 
given  fixed  nominal  plant  with  a  set  of  norm-bounded  perturbations.  If 
perturbation  is  unknown,  various  optimization  methods  (References  [80] -[83]) 
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are  applied  to  design  a  feedback  control  for  the  nominal  plant  so  that  the 
closed  loop  system  can  tolerate  "maximum"  perturbation.  If  perturbation  bound 
norms  are  given,  closed  loop  stability  conditions  are  derived  by  Doyle  and 
Stein  [69]  and  by  Chen  and  Desoer  [84]  based  on  a  generalized  version  of  the 
Nyquist  stability  theorem.  Combining  Nevanlinna-Pick  theory  and  Youla's 
parametrization,  robust  stabilizability  conditions  are  provided  by  Kimura  [85] 
for  single- input  single-output  (SISO)  systems  and  are  further  generalized  by 
Vidyasagar  and  Kimura  [86]  for  multi- input  multi-output  (MIMO)  systems. 

On  the  other  hand,  if  parameter  variation  is  of  major  concern,  then  the 
concept  of  "simultaneous  stabilization"  as  discussed  in  section  3.2  becomes 
the  appropriate  method  of  attack. 

As  mentioned  in  the  beginning  of  this  section ,  most  of  the  published 
literature  treats  these  two  major  kinds  of  uncertainties  separately.  The 
combined  presence  naturally  raises  a  question:  given  an  uncertain  system 
having  both  parameter  variation  and  unstructured  uncertainty,  is  it  possible 
to  find  a  robust  control  which  guarantees  the  closed  loop  stability  over  the 
entire  perturbation  range?  If  so,  how  to  design  a  desired  compensator?  This 
section  intends  to  examine  this  problem.  A  new  design  model  is  proposed  to 
characterize  an  uncertain  system  containing  those  two  different  uncertainties. 
Then  stability  criteria  for  closed  loop  systems  are  derived  from  classical 
results  for  this  new  model.  It  is  then  found  that  under  certain  conditions  an 
uncertain  system  represented  by  the  new  model  can  be  stabilized  by  a  proper 
(or  strictly  proper  if  desired)  and  stable  compensator.  A  procedure  of 
synthesizing  the  robust  stabilizer  is  also  provided.  All  proofs  are  in  the 
appendices . 
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Preliminary  Definitions  and  Lemmas: 

In  addition  to  the  definitions  and  lemmas  introduced  in  section  3.2,  we 
need  the  following  lemmas: 

Lemma  4 . 1  (See  Appendix  H  for  proof.) 

Assume  f(s,q),  as  in  (3.38),  is  an  n-th  order  standard  positive  Hurwitz 
invariant  uncertain  polynomial  and 

|f(jw,q)|2  A  g(w2,q)  -  ^(qjw2”  +  ^(qjw2"'2  +  ...  +  0n(q)  ,  (4.25 

then 


(1)  min(0o(q):  q  €  Q}  -  b0  >  0;  and  (4.26a) 

(2)  min(g(w2,q)  :  q  €  Q,  u  >  0)  -  k  >  0.  (4.26b) 

Lemma  4 . 2  (See  Appendix  I  for  proof.): 

Assume  f(s,q),  as  in  (3.38),  is  a  standard  positive  Hurwitz  invariant 
(PHI)  polynomial  with  degree  d(q)  and  f'(s,q)  is  a  standard  uncertain 
polynomial  with  degree  d' (q)  and  leading  coefficient  aj(q)  satisfying  the 
following  conditions:  For  all  q  €  Q, 

(1 )  d' (q)  <  d(q)  +  1,  and 

(2)  if  d'(q)  -  d(q)  +  1,  then  aj(q)  >  0. 

Then,  given  any  constant  a  >  1,  there  exists  an  c*  >  0  such  that  for  any 
c  e  (0 ,  c*) , 

fe(s,q)  -  f(s,q)  +  tf'(s,q)  (4.27a) 

is  also  a  standard  PHI  polynomial  and 


1  <  _ 

ft(jw.q)  I 

for  all  q  €  Q  and  w  >  0. 


(4.27b) 
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Modeling  of  Uncertain  Systems: 

In  the  recent  literature  on  robust  stability  theory,  various 
mathematical  models  of  uncertain  linear  systems  are  proposed  in  which, 
however,  the  two  most  widely  used  models  are  the  following: 

Model  1  (Low  Order  Parameter  Variation  Model): 

In  the  frequency  domain,  an  uncertain  system  is  represented  by  a  family 
of  rational  matrices  (P(s,q):  q  e  Q)  where  Q  is  a  prespacified  index  set,  for 
each  q  e  Q,  P(s,q)  is  an  n  x  m  rational  matrix  and  the  highest  order  of  P(s,q) 
is  a  given  finite  number. 

In  the  state  space  counterpart,  an  uncertain  system  is  represented  by  a 
family  of  matrices  (A(q) , B(q) ,C(q) ,D(q) :  q  £  Q)  in  the  state  space  equations 

x  -  A(q)x  +  B(q)U , 

(4.28) 

y  -  C(q)x  +  D(q)u,  q  €  Q, 

where  Q  is  a  prespecified  set  and  the  dimensions  of  A(q) ,  B(q) ,  C(q) ,  and  D(q) 
are  given  finite  numbers. 

Model  2  (High  Order  Norm  Bounded  Uncertainty  Model) : 

An  uncertain  system  is  represented  by  a  class  of  plants  [80] -[85] 

A(P0(s) , r(s) )  A  (P0(s)  +  L(s):  ||L(jw)||  <  |r(jw)|;  V  w  >  0)  (4.29a) 

or 

M(P0(s),r(s))  A  ((I  +  L(s))P0(s):  |L(J»)|  <  |r(jw)|;  V  «  >  0)  (4.29b) 

where  P0(s)  is  a  given  nominal  plant,  r(s)  is  a  prespecified  rational  function 
and  L(s)  is  an  unknown  but  stable  rational  matrix  whose  norm  is  bounded  by 
|r(jw)|,  A(«)  denotes  additive  model  and  M(>)  denotes  multiplicative  model. 

An  important,  but  less  recognized  fact  is  that  an  uncertain  system 
represented  by  Model  1  may  not  always  be  represented  by  Model  2  and  vice 
versa.  It  may  be  easy  to  see  that  Model  2  allows  infinite  order  uncertainty 
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Co  be  contained  while  Model  1  cannot;  it  can  also  be  seen  that  an  uncertain 
system  represented  by  Model  1  may  not  be  always  represented  by  Model  2.  lor 
example,  consider  an  uncertain  system  described  by 
P(s,q) ,  q  6  {1,2) 
in  which 

P(s,l) - ,  and 

s  -  1 


Obviously,  Model  2  is  not  suitable  for  describing  this  family  of  systems. 

Based  on  the  observations  above,  to  characterize  a  plant  containing  both 
parameter  variation  and  high  frequency  unstructured  uncertainty,  we  propose  a» 
new  model  of  uncertain  systems  combining  both  Model  1  and  Model  2  for  studying 
robust  stability  problems. 

Pef  jni  tjop_ 

An  uncertain  system  A(P0(s  ,q)  ,r(s)  ,Q)  is  defined  by 
A(P0(s,q),r(s),Q)  A  {P0(s,q)  +  L(s):  q  6  Q;  ||L(j«)fl  <  | r(jw) |  ,  V  w  >  0}  (4.30) 

and  an  uncertain  system  M(P0(s,q) ,r(s) ,Q)  is  defined  by 

M(P0(s,q),r(s),Q)  A  {(I  +  L(s)  )P0(s  ,q)  :q  e  Q;flL(jw)||  <  |r(jw)  |,  V  «  >  0)  (4.31) 

where  Q  is  a  prespecified  index  set;  for  each  q  £  Q,  the  nominal  plant  P0(s,q) 
is  an  n  x  m  rational  uncertain  matrix;  r(s)  is  a  prespecified  rational 
function  and  L(s)  is  any  unknown  but  stable  rational  matrix  whose  norm  is 
bounded  by  |r(jw)  |  . 

Stability  Criteria: 

In  the  previous  section,  we  proposed  a  new  model  of  uncertain  systems. 

Our  main  objective  is  to  design  a  robust  feedback  control  to  stabilize  an 
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uncertain  system  represented  by  the  new  model.  To  this  end,  we  need  stability 
criteria  for  the  closed  loop  systems.  First,  we  introduce  a  standard  result 
for  uncertain  systems  with  only  unstructured  uncertainty  [86], 


Lemma  4 . 3 

If  a  controller  C(s)  stabilizes  P0(s) ,  then  C(s)  stabilizes  all  plants 
in  the  class  A(P0(s),  r(s))  if  and  only  if 
| C ( s ) ( I  +  P0(s)C(s) )'1r(s) I,  <  1. 

This  lemma  can  be  naturally  extended  for  our  new  model;  no  proof  is  needed. 
Theorem  4. 1 

If  C(s)  stabilizes  P0(s,q),  then  C(s)  stabilizes  all  plants  in  the  class 
A(P0(s,q) ,r(s) ,Q)  if  and  only  if 

II C ( s )  ( I  +  P0(s,q)C(s))-1r(s)||-  <  1 
for  all  q  6  Q. 

For  SISO  systems,  a  straightforward  computation  directly  leads  to  the 
following  corollary. 

SsxisILary.AJi 

A  controller  C(s)  A  Ne(s)/Dc(s)  stabilizes  a  SISO  uncertain  system 
A(P0(s,q) ,r(s) ,Q)  with  P0(s,q)  A  N(s,q)/D(s,q)  if 

A(s,q)  A  De(s)D(s,q)  +  Nc(s)N(s,q)  (4.32) 

is  Hurwitz  invariant  over  Q  and 

|N<i(Jw)N(Ja>,q)|  |r(Jw)|  (4  33) 

|Dc(j«)D(j«,q)  +  Ne(j«)H(J«,q)  |  “ 

for  all  q  €  Q  and  u  >  0. 

Robust  Stabllizabllitv  of  SISO  Uncertain  Systems: 

Consider  a  SISO  uncertain  system  A(P0(s,q) ,r(s) ,Q)  or  M(P0(s,q) ,r(s) ,Q) 
as  in  Definition  4.1,  where 
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m(q) 

l  °i(q)s®(<,>  i 

P  (s  q)  A  a  _ 

o(  ’q;  *  D(s,q)  A  m'(q) 

l  «;(q)«",<,,_1 

i-0 


where  the  symbols  have  the  same  definitions  as  in  (3.40).  The  following 
assumptions  play  a  central  role  for  robust  stabilizability . 

Assumption  4.1  (Standardness) 

P0(s,q)  is  a  standard  uncertain  rational  function. 

Assumption  4.2  (Minimum  Phase) 

For  each  fixed  q  €  Q,  P0(s,q)  is  of  minimum  phase;  i.e.,  all  zeros  of 
P0(s,q)  lie  in  the  strict  left  half  plane. 

Assumption  4.3  (One  Sign  High  Frequency  Gain) 

The  sign  of  the  ratio  of  the  leading  coefficients  of  D(s,q)  and  N(s,q) 
is  sign  invariant  over  Q.  In  the  sequel,  we  keep  the  sign  of  the  leading 
coefficient  of  D(s,q)  positive. 

Assumption  4.4  (Boundedness) 

For  an  uncertain  system  A(P0(s,q) ,r(s) ,Q) , 

SUP  (iwT  qeQ;  »>oJ<l. 

For  an  uncertain  system  M(?0(s ,q) , r(s) ,Q) , 
sup{  |r(jw)  | :  u>  >  0)  <  1. 

Remark:  An  uncertain  system  satisfying  Assumptions  4. 1-4.4  can  be 
viewed  as  a  family  of  finite  or  infinite  number  of  plants.  Each  plant  may 
have  a  different  nominal  plant  with  a  different  relative  degree  and  unknown 
but  norm  bounded  stable  high  frequency  uncertainty.  The  zeros  of  each  nominal 
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plant  are  in  the  strict  left  half  plane.  There  is  no  restriction  on  pole 
positions  of  the  nominal  plants. 

Theorem  4 , 2  (See  Appendix  J  for  proof.) 

Under  the  Assumptions  4. 1-4.4,  a  SISO  uncertain  system  A(P0(s ,q) , r(s) ,Q) 
or  M(P0(s ,q) ,r(s) ,Q)  is  robust  stabilizable.  Furthermore,  a  robust  stabilizer 
C(s)  can  be  constructed  to  be  proper  (or  strictly  proper  if  desired)  and 
stable.  The  minimum  order  of  C(s)  is  equal  to 

k  -  max(m'(q)  -  m(q) :  q  e  Q)  -  1  or  k  -  0  when  max(m* (q)  -  m(q) :  q  e  Q)  -  0. 

Remark :  When  high  frequency  unstructured  uncertainty  is  neglected,  that 
is,  L(s)  ■  0,  Theorem  4.2  coincides  with  the  main  result  in  Reference  [60]. 
When  parameter  variations  are  not  under  consideration,  i.e.,  P0(s,q)  reduces 
to  PQ(s) ,  it  is  of  interest  to  compare  Theorem  4.2  with  some  existing  results, 
such  as  References  [69]  and  [85].  First,  both  references  require  a  full  order 
controller,  while  Theorem  4.2  provides  a  significantly  lower  order 
compensator,  depending  only  on  the  largest  relative  degree  of  all  nominal 
plants.  Secondly,  the  stability  margin  for  a  system  M(P0  (s) ,r(s) ,Q)  with  a 
minimum  phase  nominal  plant  is  improved  to  inf|r(jw)|  <  1  from  |r(jw)  |  <  1/2 
obtained  by  applying  LQG/LTR  in  Reference  [69].  Finally,  in  Theorem  4.2  there 
are  no  special  restrictions  on  the  relative  degree  of  r(s)  and  the  locations 
of  unstable  poles  of  P0(s) ,  as  in  Reference  [85]. 

The  Stabilizing  Compensator  Synthesis: 

In  this  section,  we  extract  the  crucial  ingredients  for  compensator 
construction  from  the  proof  of  Theorem  4.2. 

Procedure  of  Synthesis  4.1 


Step  1:  Consider  a  compensator  of  the  form 


Ne(s) 


Avs*  +  A, 


„k-l 


+  ...  +  A,s  4-  An 


C(s)  A  A  ,  , 

Dc(s)  CjS1  +  Cj^s1  1  +  ...  +  C1s  +  £0 


1  >  k. 


•Is  '  ‘i-ia 

Choose  Nc(s)  to  be  any  Hurwitz  polynomial  with  Ak  >  0  when  N(s,q)  is  PHI  or 
Ak  <  0  when  N(s,q)  is  negative  Hurwitz  invariant  and  degree 
k  >  max{m'(q)  -  m(q):  q  6  Q)  -1 
or  k  -  0  when  max{m'(q)  -  m(q) :  q  €  Q)  -  0. 

Step  2  (Initialization):  Define 
A0(s,q)  A  Nc(s)N(s,q) ; 

Dc,0(s)  A  0; 

IPoO-q) 


(4.34) 


5  -  inf 


q  6  Q;  “  >  o  [  ■  i 


|r(»| 

for  an  uncertain  system  A(P„(s , q) , r(s) ,Q)  and 


6  -  inf  i 


U<»  I  : 


q  e  Q; 


>  u  >  0  |  * 


for  an  uncertain  system  M(P0(s,q) ,r(s) ,Q) . 

S£gP  3  (Inductive  Step):  Given  ^(s.q)  and  Dc ^(s),  select  ct  >  0  such 

that 


A1+l(s.q)  A  A1(s  ,q)  +  c^DCs.q)  (4.35) 

is  Hurwitz  invariant, 

Dc  l+1(s)  A  DCii(s)  +  (4.36) 

is  Hurwitz  (when  a  stable  controller  is  required)  and 

|At(j«,q)  I  1  +  S/2  +  S/k  +  ...  +  S/2 in 
Ui+i(jw,q)  |  "  1  +  S/2  +  5/4  +  ...  +  5/21  (4‘37) 

for  all  q  €  Q  and  w  >  0.  Assumptions  4. 1-4.4  guarantee  the  success  of  this 

procedure  for  i  -  0,1,2 . 1. 

Step  4  (Termination) :  The  denominator  of  the  compensator  is  given  by 
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Dc(s)  A  Dc<1+1(s). 

Remark:  With  the  model  proposed  in  the  last  section,  one  may  think  of 
an  alternative  design  procedure  for  a  robust  stabilizer:  take  a  plant  P„(s) 
having  the  highest  relative  degree  from  P0(s,q)  as  the  nominal  plant.  Design 
a  controller  to  stabilize  P0(s)  with  its  high  frequency  uncertainty  and  then 
increase  the  control  gain  to  stabilize  the  rest  of  the  systems  in  the  family. 
The  following  example  shows  that  this  method  cannot  lead  to  a  desired  robust 
stabilizer,  even  for  an  uncertain  system  without  high  frequency  unstructured 
uncertainty.  For  example,  consider  a  family  of  systems: 


100 

P0(s,l)  -  -Jr; 

_  ,  s2  +  2s  +  1 

Vs -2)  “  s*  .  20s3  • 

Obviously,  the  compensator 

r,  ^  s2  +  2s  +  1 

(S)  s2  +  10s  +  50 

stabilizes  the  system  P(s,l).  However,  it  is  easy  to  verify  that  for  any 
k  €  R,  the  controller  C1(s)  -  kC(s)  cannot  stabilize  the  system  P(s,2). 
However,  with  the  method  given  in  [60],  one  may  find  a  robust  stabilizer  as 
follows : 


.  2  x  1q3(s2  +  2s  +  1) 

(S  “  s2  +  8000s  +  50000 

Illustrative  Examples 

In  this  section,  two  examples  are  provided  to  illustrate  the 
construction  of  a  robust  compensator  for  the  combined  uncertainty  case. 
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Example  4.1 

An  uncertain  system  M(P0(s,q) ,r(s) ,Q) ,  as  in  (4.31),  is  characterized  by 
Q  -  {1,2,3}  and 


P0(s,l) 

P0(s,2) 

P0(s,3) 


1 

s2  -  2s  -  2’ 
1 

s  -  1  ' 

s2  +  2s  +  1 
s3  -  2s2  : 


r(s) 


0.1 


s  +  1 

It  is  easy  to  verify  that  Assumptions  4. 1-4.4  are  satisfied;  hence,  the  system 
is  stabilizable.  A  robust  compensator  C(s)  can  be  constructed  as  follows: 

Step  1:  We  now  require  k  >  max{2,l,l)  -  1  -  1.  Taking  k  -  1  and  1-1, 
the  compensator  has  the  form 
A,s  +  A. 

C(s) - 1  — A,  >  0. 

£is  +  co 

Next,  we  choose  Axs  +  A0  -  s  +  1,  i.e.,  Ne(s)  -  s  +  1. 

Step  2 :  We  have 
A0(s,l)  -  s  +  1, 

Vs, 2)  -  s  +  1, 

A„(s,3)  -  (s  +  l)(s2  +  2s  +  1),  and 
&  -  inf  { 1/ | r(jw)  | :  w  >  0)  -  l  -  9. 

Steps  3-4:  To  select  e„,  consider 
Aj(s , 1)  -  £0(s2  -  2s  -  2)  +  Afl ( s , 1 ) , 

A1(s  ,  2)  -  £0(s  -  1)  +  A0 ( s ,  2 ) , 

A1(s , 3)  -  £0(s3  -  2s2)  +  A0(s,3) 

and 
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KO,q)  | 

Ui(J«.q)  I 


<1+5/2 


5.5  for  q  6  {1,2,3} . 


In  view  of  all  three  cases,  we  can  take  c„  -  0.1  to  assure  Hurwitz  invariance 
of  A1(s,q)  and  satisfaction  of  the  inequality  above.  Now,  we  again  consider 
three  cases  of 

A2(s,l)  -  s(s2  -  2s  -  2)  +  A^s.l), 

A2(s,2)  -  s  -  1)  +  A1(s , 2) , 

A2(s,3)  -  £xs(s3  -  2s2)  +  Al(s,3), 

and 


KCKq)  I  1  +  5/2  +  5/4 
U2(j«,q)  |  -  1  +  5/2 


1.409  for  all  q  e  (1,2,3). 


we  can  assure  Hurwitz  invariance  of  A2(s,q)  and  satisfaction  of  the  inequality 


above  by  taking  £x  -  0.005.  Hence,  a  stabilizing  compensator  is  given  by 

s  +  1  200(s  +  1) 

(S)  ”  0.005s  +  0.1  “  s  +  20  ‘ 


Example  4,2 


An  uncertain  system  M(P„(s ,q) ,r(s) ,Q) ,  as  in  (4.34),  is  characterized  by 


p0<s.q) 


_  Si. 
s  +  q2 


ql  6  [1,2] ;  q2  e  [-1,1] , 


In  this  example,  the  uncertain  system  has  continuous  parameter 
variation.  It  is  easy  to  verify  that  Assumptions  4. 1-4.4  are  satisfied, 
hence,  the  system  is  stabilizable.  A  robust  compensator  C(s)  can  be 
constructed  as  follows: 

Step  1:  We  now  require  k  >  0.  Taking  k  -  0  and  1-1,  the  compensator 
has  the  form 


79 


C(s)  - 


Next,  we  choose  X0  -  1,  i.e.,  Ne(s)  -  1. 
Step  2 ;  We  have 

Vs.q)  -  qi 

6  -  Inf (1/ | r(j«)  |  :  w  >  0)  -  1  -  1. 


Steps  3-4:  Now 


A^s.q)  -  Cq( s  +  q2)  +  -  e„s  +  (c0q2  +  qt)  . 


Ko.q)  I 

h1(j<*>.  q) 


<  +6/2; 


or 

+  (q:  +  c0q2)2 "  9/4 

With  the  given  bounds  of  q4,  we  require  c0  <  1  for  Hurvitz  invariance  of 
A1(s,q)  and  c0  <  1/3  for  satisfying  the  inequality  above.  Hence,  we  take 
e0  -  0.2.  A  straightforward  computation  results  in 
A2(s,q)  -  c1s2  +  e^s  +  0.2s  +  0.2q2  +  qx 


and  we  require 


iMJw.q)  I  <  1  +  S/2  +  S/4 

MjwTq)  |  ~  l  +  s/2 


or 


0.04^2  +  (q,  +  0.2q,)2 _ 

(cxq2  +  0.2)2  +  (0.2q2  +  qt  -  c^2)2  ~  / 

With  the  given  bounds  of  qt,  we  can  assure  Hurvitz  invariance  of  A2(s,q)  with 

ct  <  0.2  and  satisfaction  of  the  inequality  above  with  cx  <  0.0026.  Hence,  we 

can  take  c:  -  0.002  and  the  stabilizing  compensator  is  now  given  by 

1  500 


C(s)  - 


0.002s  +0.2  s  +  100  ' 
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In  this  section,  we  proposed  a  new  design  model  for  uncertain  systems 
containing  both  parameter  variation  and  unstructured  uncertainty.  Sufficient 
conditions  are  provided  under  which  robust  stabilizability  of  SISO  uncertain 
systems  is  guaranteed.  Furthermore,  the  robust  compensator  can  be  constructed 
to  be  strictly  proper  and  stable;  its  numerator  can  be  any  stable  polynomial 
having  degree  k  given  in  Step  1  of  Procedure  4.1.  Its  denominator  has 
coefficients  which  can  be  computed  one  at  a  time  using  the  recursion  given 
in  Step  3  of  Procedure  4.1. 

Some  Comments  on  Computational  Requirements: 

The  proposed  design  procedure,  of  course,  demands  a  fairly  good  amount 
of  computation,  both  in  terms  of  labor  and  complexity.  Firstly,  one  needs  an,, 
algorithm  to  test  the  Hurwitz  invariance  of  uncertain  polynomials.  The 
Kharitonov  test  and  related  improvements  as  discussed  in  section  H  would  be 
of  assistance  in  this  regard.  Next,  one  needs  to  evaluate  the  inequality  of 
(4.37)  for  all  frequencies  w  >  0.  This  step  is  more  laborious  than  complex. 
One  attractive  or  comforting  feature  of  this  design  procedure  is  that  it  is 
recursive  in  nature.  Thus  at  any  given  step  in  the  computational  procedure, 
there  is  only  one  parameter  to  be  determined,  which  significantly  simplifies 
the  computational  burden. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH  ] 

i 

I 

5.1  THE  WORK  IN  RETROSPECT 

The  main  theme  of  the  described  research  under  the  present  contract  has 
been  to  analyze  and  synthesize  controllers  for  robust  stability  for  linear 
time  invariant  systems  with  structured  parametric  uncertainty,  both  from  a 
time  domain  and  a  frequency  domain  perspective.  First  the  aspect  of  analysis 
for  stability  robustness  is  considered.  A  useful  review  of  currently 
available  methods  of  analysis  (namely,  "polynomial"  and  "matrix"  methods)  is 
given.  Then,  an  algorithm  is  presented  to  reduce  the  conservatism  of 
Kharitonov  test  for  testing  the  Hurwitz  invariance  of  uncertain  polynomials 
with  dependent  coefficients.  Also,  for  systems  described  by  state  space 
models,  improved  upper  bounds  for  robust  stability  are  obtained  by  further 
exploiting  the  functional  dependence  of  the  elements  of  the  perturbation  on 
the  uncertain  parameters.  The  contents  of  this  section  cover  the  task 
described  as  "task  3"  in  the  original  proposal. 

Next,  the  aspect  of  synthesis  of  controllers  is  addressed.  A  robust 
control  design  method  is  presented  for  systems  described  by  state  space 
representation  such  that  it  maximizes  (in  some  sense)  the  stability  robustness 
bound  for  a  given  structure  of  the  uncertainty,  taking  into  account  the 
control  effort  constraints.  This  algorithm  adequately  answers  the  question 
raised  in  "task  1"  of  the  proposal.  Then,  in  addition,  another  design  method 
from  frequency  domain  point  of  view  is  also  presented  such  that  a  single 
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compensator  simultaneously  stabilizes  a  class  of  non-minimum  phase  systems. 

The  conditions  under  which  the  existence  of  the  controller  is  guaranteed  are 
given  along  with  an  algorithm  to  construct  the  compensator. 

Then  the  aspect  of  combined  structured  and  unstructured  uncertainty  is 
considered.  First  in  the  time  domain  framework,  specifically  for  large  space 
structure  models,  a  problem  formulation  is  presented  in  which  the  unmodeled 
dynamics  are  characterized  as  an  additive  perturbation  and  control  gain  bounds 
are  derived  such  that  the  closed  loop  system  remains  stable  even  in  the 
presence  of  both  control  and  observation  spillovers.  The  design  algorithm  is 
then  extended  to  account  for  parameter  variations  in  the  lower  order  control 
design  model.  In  other  words,  the  unstructured  and  structured  uncertainty 
cases  were  obtained  in  sequence  of  one  another  rather  than  being 
simultaneously  present.  However,  the  design  algorithm  can  be  extended  to  the 
case  when  both  types  of  errors  are  present  simultaneously.  This  section  of 
the  report  covers  the  task  defined  as  "task  2"  of  the  original  proposal. 

Finally,  the  aspect  of  simultaneous  presence  of  structured  and 
unstructured  uncertainty  in  the  system  is  addressed  from  frequency  domain 
viewpoint.  An  algorithm  is  presented  to  design  a  single  (reduced  order) 
compensator  that  guarantees  robust  stabilization  under  the  combined  presence 
of  these  two  types  of  uncertainties .  The  conditions  for  the  existence  of  the 
controller  and  its  construction  are  also  given  in  detail.  This  result  covers 
the  task  raised  as  "task  4”  of  the  proposal. 

The  publications  listed  as  References  [87] -[97]  are  the  result  of  this 
research  study. 
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As  it  normally  occurs,  another  result  of  this  study  is  that  many 
interesting  research  topics  have  surfaced  for  further  investigation.  These 
are  described  in  the  next  section. 

5.2  TOPICS  FOR  FURTHER  RESEARCH  WHICH  NEED  CONTINUED  SUPPORT 

(  1)  The  foremost  area  of  research  would  be  to  extend  the  frequency  domain 
design  methodology  for  the  combined  structured  and  unstructured  uncertainty 
problem  by  relaxing  some  of  the  assumptions  made  with  respect  to  the  minimum 
phase  condition,  and  to  enlarge  the  class  of  unstructured  uncertainty  profiles 
considered  in  the  synthesis  of  the  compensator  for  robust  stabilization. 

(2)  Another  topic  for  further  research  would  be  the  incorporation,  in 

ft- 

addition  to  the  robust  stabilization  requirement,  of  the  requirement  of 
performance  robustness  in  the  problem  formulation.  Some  specific  constraints 
to  be  considered  could  be  disturbance  rejection  and  bounded  control  effort. 

(  3)  An  interesting  area  of  research  would  be  to  develop  the  state  space 
conterpart  to  the  frequency  domain  design  methodology  given  for  the  combined 
uncertainty  problem. 

(  4)  One  aspect  that  needs  attention  is  the  extension  of  the  proposed 
frequency  domain  methodology  under  combined  uncertainty  problem  to  Multiple 
Input  Multiple  Output  (MIMO)  systems.  This  is  definitely  not  a 
straightforward  extension  of  the  Single  Input  Single  Output  (SISO)  case. 

(  5)  It  is  also  important  to  explore  the  possible  relationship  of 
different  methodologies  available,  such  as  structured  singular  value  design, 
multivariable  stability  margin  design,  simultaneous  stabilization  design  using 
polynomial  theory,  and  the  H*  control  design,  since  all  of  these  approaches 
seem  to  solve  the  same  problem. 
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(  6)  Of  course,  it  is  always  beneficial  to  apply  the  developed 
methodologies  to  specific  applications,  such  as  aircraft/spacecraft  control, 
control  of  flexible  structures,  and  robotics. 
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APPENDIX  A:  PROOF  OF  THEOREM  2.1 


Necessity  is  trivially  proved  by  letting  kx(q)  -  k^q)  “  1. 

(Sufficiency) :  When  n  is  even,  let 

h(s2,q)  -  a0(q)sn  +  a2(q)sn‘2  +  o*(q)sB'*  +  ...  +  aB(q) 

and 

sg(s2,q)  -  a1(q)s“*1  +  o3(q)sn‘3  +  a5(q)sn*5  +  ...  +  a^.^qjs. 

Then 

f'(s.q)  -  k1(q)h(s2,q)  +  k^sg^.q) 

-  k1(q)h(z,q)  +  sk2(q)g(z,q). 

It  follows  from  the  Hermite-Bieler  theorem  that  if  f'(s,q)  is  Hurwitz 
invariant  over  Q,  then  k^qWz.q)  and  kj(q)g(z,q)  form  a  positive  pair  for 
every  q  €  Q,  i.e.,  the  zeros  of  the  polynomials  k1(q)h(z,q)  and  k^q^z.q) 
must  be  distinct,  real,  negative,  and  interlaced  as  follows: 
ua  <  v2  <  Uj  <  v2  <  ...  <  <  u.  <  0 

where  m  —  n/2,  uA  are  zeros  of  k2(q)h(z,q)  and  vt  are  zeros  of  k2(q)g(z,q). 
However,  kx(q)h(z,q)  and  h(z,q)  have  the  same  zeros.  Hence,  h(z,q)  and  g(z,q) 
also  form  a  positive  pair  for  every  q  €  Q.  Thus,  the  Hermite-Bieler  theorem 
guarantees  that 

f(s,q)  -  h(z , q)  +  »g(z,q) 
is  also  Hurwitz  invariant  over  Q. 

When  n  is  odd,  a  similar  proof  leads  to  the  same  conclusion. 
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APPENDIX  B:  PROOF  OF  THEOREM  2.3 


By  assumption  the  four  bounding  polynomials  of  f(s,q)  are  Hurwitz.  It 
follows  from  the  Kharitonov  theorem  that  every  polynomial 
f(s)  -  a0sn  +  atjs”'1  +  ...  +  an,  aL  €  [a1,blJ 
is  Hurwitz.  Since 

ai  <  ai-  bI  <  bi> 

then  every  polynomial 

f(s)  -  a0sn  +  c^s"*1  +  ...  +  an>  al  €  [a^.bj] 
is  also  Hurwitz.  Consequently,  the  four  bounding  polynomial?  of  f'(s,q)  are 
Hurwitz . 
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APPENDIX  C:  PROOF  OF  THEOREM  2.4 


IC  is  known  Chat  Che  perturbed  system 
x(t)  -  (A0  +  E)x(t) , 

where  A0  is  an  asymptotically  stable  matrix,  is  asymptotically  stable  if 

sup  p((jwl  -  Ag)‘lE]  <  1 
o)>0 

Let  (jwl  -  A0)_1  -  M(w) .  Now  with  E  given  by 
r 

E  -  £  fit Et, 
i-1 

the  perturbed  system  (F.l)  is  asymptotically  stable  if 

sup  p  [  M(w)  [  X  <  1 

«>0  L  M-l 

But 

max  |^j|  sup  p  [  X  |M(«)El  ]  >  sup  p  [  X  I^MO^Ej  1 

j  w>0  i-1  w>0  '•i-1  J 

>  ®up  p  [  X  /9iM(«)E1  1 

w>0  '■i-1  J 

r 

>  sup  p  [  M(w)(  X  P^i)  ] 

w>0  i-1 

The  satisfaction  of  condition  (2.12a)  implies  the  satisfaction  of  (C.4). 

Hence,  the  perturbed  system  is  asymptotically  stable. 

For  r  -  1,  we  see  that 

sup  p  [M(w)E1)J  >  sup  p  [M(w)  (/JjEj)  ] 

w>0  w>0 

Hence  (2.12b)  implies  (C.4),  leading  to  the  result. 
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APPENDIX  D:  PROOF  OF  LEMMAS  3.3  AND  3.4 


By  the  hypothesis  of  Lemma  3.3,  the  index  set  Q  can  be  divided  into  a 
finite  number  of  subsets  : 

Qmr  ■  (q  e  Q:  d(q)  “  n;  d' (q)  -  n' ) . 

In  each  subset  Q^, ,  it  follows  from  Lemma  3.1  that  there  exists  an  c^.  >  0 
such  that  for  any  c  €  (O.t^,), 

f£(s .q)  -  f(s,q)  +  £ f * (s ,q) 

is  also  PHI  over  Q^, .  Hence,  the  conclusion  of  Lemma  3.3  immediately  follows 
by  taking 


A  similar  proof  applies  to  Lemma  3.4. 
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APPENDIX  E:  PROOF  OF  THEOREM  3.2 

By  assumption  that  N^s.q)  is  Hurwitz  invariant,  it  follows  from  Theorem 
3.1  that  there  exists  a  C0(s)  -  Ne(s)/De  „(s)  such  that 
A„(s,q)  -  N^sJN^s.q)  +  De  0(s)D(s,q) 
is  PHI  over  Q. 

Now,  applying  Lemma  3.4,  there  exists  an  c0  >  0  such  that 
A^s.q)  -  sA0(s , q)  +  c0De  0(s)D(s ,q) 
is  PHI  over  Q.  Similarly,  there  exists  an  ct  >  0  such  that 
A2(s,q)  -  sA1(s , q)  +  c1De>0(s)D(s ,q) 
is  PHI  over  Q,  and  there  exists  an  c2  >  0  such  that 
A3(s,q)  -  sA2(s,q)  +  c2De  0(s)D<S,q) 

-  s3A0<s.q)  +  (c0sa  +  £xs  +  ea)Dc  0(s)D(s,q) 
is  PHI  over  Q  and  Dc  3(s)  -  £„s2  +  £:s  +  e2  is  Hurwitz. 

In  the  next  step,  it  follows  from  Lemma  3.4  that  there  exists  an  c*  >  0 
such  that  for  any  c  €  (0,c*), 

fc(s,q)  -  sA,(s)  +  De  „(s)D(s,q) 

is  PHI  over  Q.  Also,  there  exists  an  c**  >  0  such  that  for  every  c  e  (0,e**), 
fj(s)  -  sDc  3(s)  +  t 

is  Hurwitz.  Hence,  we  can  take  0  <  c3  min(c*,c**)  to  guarantee  that 
A*(s,q)  -  sA3(s)  +  e3Dc  0(s)D(s,q) 
is  PHI  and 

De  4(S)  -  SDC  3(S)  +  C3  -  £„S3  £ XS2  +  £jS  +  £j 

is  Hurwitz. 

Repeating  the  above  procedure  for  X*1  times,  we  complete  the  proof. 
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APPENDIX  F:  PROOF  OF  THEOREM  3.3 

By  applying  Lemma  3.4,  we  take  f(s,q)  -  Dx(s,q)  and  f'(s,q)  -  N(s,q)  and 
follow  the  same  procedure  as  in  Appendix  E  to  construct 
Ne(s)  -  60sx'1  +  S^'2  +  ...  +  Sw 
such  that 

Ai(s,q)  -  Nc(s)N(s,q)  +  s^s.q) 
is  PHI.  Then  we  apply  Lemma  3.3  by  taking  f(s,q)  -  A£(s,q)  and 
f ' (s , q)  -  sD(s,q)  -  sA+1D1(s,q).  Hence,  there  exists  an  tj  >  0  such  that 
&'2(s,  q)  -  AJ(s,q)  +  ^sDCs.q)  -  Ne(s)N(s,q)  +  (1  +  s)D(s,q) 
is  PHI  over  Q.  Similarly,  there  exists  an  c2  >  0  such  that 
A'(s,q)  -  A'(s,q)  +  <2s2D(s,q) 
is  PHI  over  Q  and  Dc  2(s)  -  1  +  CjS  +  e2s2  is  Hurwitz. 

In  the  next  step,  it  follows  from  Leona  3.4  that  there  exists  an  t*  >  0 
such  that  for  any  e  €  (0,c*\ 

ft(s,q)  -  Aj (s)  +  cs3D(s,q) 

is  PHI  over  Q.  Also  there  exists  an  t **  >  0  such  that  for  every  c  e  (0,c**), 
f;(s)  -  De  2(s)  +  cs3 

is  Hurwitz.  Hence,  we  can  take  0  <  c3  <  nin(e*,e**)  to  guarantee  that 
A;(a,q)  -  AJ(s)  +  £3s3D(s,q) 
is  PHI  and 

De  3(s)  -  Dc  2(s)  +  c3s3  -  1  +  c:s  +  e2s2  +  e3s3 
is  Hurwitz. 

Repeating  the  above  procedure  for  1  >  A- 1  tines,  we  complete  the  proof. 
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APPENDIX  G:  PROOF  OF  THEOREM  3.4 

By  assumption  that  P^s.q)  is  simultaneously  stabilizable ,  there  exists 

a  compensator  C^s)  -  No  1(s)/Dc  ^s)  such  that 

A(s,q)  -  Dcl(s)D1(s,q)  +  Ne  1<s)N1(s ,q) 

is  Hurwitz  invariant  over  Q.  Now,  construct  a  compensator 

c(s)  . . . -  .  ■  Nc.i(«)Dc(^> 

1  -  C1(s)P0(s)  DQ(s)De  l(s)  -  N0(s)Ne  l(s) 

and  the  closed  loop  system 

.  P(s.q)  N(s,q)Dr(s) _ 

,q;  “  1  +  C(s)P(s,q)  "  D(s,q)Dc(s)  +  N(s,q)Nc(s) 

_ N(s,q)Df (s) _ 

Dg(s)D1(s,q)De  ^(s)  +  N^s.qJN,  ^(s) 

Since  D0(s)  is  Hurwitz  and  D1(s,q)De|1(s)  +  Nj(s , q)Nc  j(s)  is  Kurwitz  invariant 
over  Q,  the  closed  loop  system  is  robust  stable. 
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APPENDIX  H:  PROOF  OF  LEMMA  4.1 

Knowing  that  j90(q)  -  a2(q)  is  a  continuous  function  in  the  compact  set  Q 
and  a0(q)  >  0  for  all  q  G  Q,  it  immediately  follows  that 
min  (0o(q) :  q  6  Q)  -  b„  >  0. 

Since  f(s,q)  is  PHI,  we  have 
g(w2,q)  >  0 

for  every  q  6  Q  and  u>  >  0.  The  standardness  of  f(s,q)  guarantees  that 

0o<q).0i<q> . 0„<q)  are  a11  bounded  and  min  {0o(q):q  €  Q)  -  b0  >  0,  so  for 

any  kx  >  0,  there  exists  an  (^  >  0  such  that  for  all  u  >  and  q  6  Q, 
g(w2,q)  >  kx. 

Now  g(w2,q)  is  a  continuous  function  in  a  compact  set  so  it 

follows  that 

min  (g(w2,q):  q  e  Q,  w  6  [O.wJ)  -  I4  >  0. 

Consequently , 

min  (g(w2,q):  q  €  Q,  u  >  0}  -  min  (kj.kj)  -  k  >  0. 
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APPENDIX  I:  PROOF  OF  LEMMA  4.2 

By  assumption  that  f(s,q)  and  f'(s,q)  are  standard,  the  index  set  Q  can 
be  divided  into  a  finite  number  of  compact  subsets  Q^, ,  where 
Qnn*  -  (q  €  Q:  d(q)  -  n,  d' (q)  -  n' } . 

First,  we  prove  that  the  lemma  is  true  when  q  €  i.e.,  given  any  a  >  1, 

there  exists  an  e^,  >  0  such  that  for  any  e  e  (0,^,), 
fe(s,q)  -  f (s , q)  +  cf'(s,q) 
is  also  PHI  over  and 

I 

for  all  q  €  Q^,  and  w  >  0.  Then,  by  letting 
c*  -  min  (e^.:  V  n.n' } , 
the  conclusion  of  Lemma  4.2  directly  follows. 

It  follows  from  Lemma  4.5  that  there  exists  an  cx  >  0  such  that  for  any 
c  6  (0,cx) , 

fe(s,q)  -  f(s,q)  +  ef’(s,q) 
is  also  PHI  over  Q^, . 

We  now  claim  that  for  any  a  >  1,  there  exists  an  c2  >  0  such  that  for 
any  c  €  (0,c2), 

^J.<. 

«.q)  I 

for  all  q  €  Q^,  and  w  >  0.  Once  the  claim  is  proved,  the  conclusion 
immediately  follows  by  letting  cm,  -  min  {c1,e2}. 


Ill 


f 


To  prove  the  claim,  it  is  equivalent  to  show  that  there  exists  an  cz  >  0 
such  that  for  any  c  e  (0,c2), 

<  1  +  6X 

for  all  q  e  and  w  >  0  where  Sx  -  a2  -  1  >  0;  or,  equivalently, 

(1  +  Sx)  | fe(Jw,q)  |2  -  I f (Jw.q)  |2  >  0. 

Case  1:  n'  -  n  +  1 

A  straightforward  computation  yields 
(1  +  Sx)  |fe(J«,q)  |2  -  | f (jw,q)  j2 

-  (1  +  $1)c2aJ2(q)w2a+2  +  hjtw2^, +  «x|f(Jw,q)|2 

-  (1  +  fi^c^^CqJu2"*2  +  hx(»2,q,c,«x)  +  «xg(w2,q) 
d  (1  +  «1)c2aJ2(q)w2n+2  +  h(«2,q,e,*x) 

where  aj2(q)  is  the  leading  coefficient  of  f'(s,q)  and 

h1(w2,q,e,fil)  -  ciygCq.c.SJu2*  +  yxtq.c.fiiJw2”'2  +  ...  +  7n(q, £,«!>) 
g(w2.q)  A  |f(jw,q)  |2  -  ^0(q)«2B  +  ^1(q)w2"’2  +  ...  +  0„(q)  . 

The  standardness  of  f(s,q)  and  f'(s,q)  guarantees  that  for  given  $x  when 
c  ■*  0,  c71(q,c,«1)  uniformly  converges  to  zero  in  the  compace  set  for  all 
i  -  0,1,..., n.  By  Lemma  4.1, 

min  I0o(q):  q  €  Q^, )  -  b0  >  0. 

So,  there  exists  an  e3  >  0  such  that  for  any  c  €  (0,c3), 
min  {«x0o(q)  +  C70(q,e,«x):  q  €  Q^,.)  >  0.5«xb0  >  0. 

Consequently,  for  any  Mx  >  0,  there  exists  a  ua  >  0  such  that  for  any  u  >  ua 


1  f  ( j<*>.  q)  1 2 

I2 


and  q  €  Q^, , 

h(w2,q,c,«1)  >  Mx. 
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For  given  Slt  h(«2,q,c,51)  is  s  continuous  function  of  c  in  the  compact  set 
{ [0.wB]  .Qnn* )  >  it  follows  that  when  t  -*  0,  h(w2,q,c,il)  uniformly  converges  to 
51g(w2,q).  Hence,  there  exists  an  e4  >  0  such  that  for  any  e  6  (0,c4), 
min  (h(w2,q,£,51) :  q  e  Q^,, ,  w  6  [O.uJ)  >  0.5«xk  >  0 
where  k  -  min  (g(w2,q):  q  e  «  >  0)  >  0  by  Lemma  4.6.  Now 

take  c2  -  min  {c3,£4)  >  0  and  M  -  min  (M1,0.562k}.  Then,  for  any  e  e  (0,c2), 
min  (h(£2,q,  c  ,6^  :  q  €  Q^. ,  w  >  0}  >  M  >  0. 

Since  (1  +  61)£2aj2(q)w2n^2  >  0,  we  conclude  that  for  any  c  €  (0,c2), 

(1  +  f^If.O.q)!*  -  |f(j«,q>|2  >  0 
for  all  q  e  Q^,.  and  w  >  0. 

Case  2:  n'  <  n  +  1 

In  this  case,  the  term  w2”*2  no  longer  exists  and  the  proof  above  remains 

valid. 
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APPENDIX  J:  PROOF  OF  THEOREM  4.2 


Consider  an  uncertain  system  A(P0(s ,q) ,r(s) ,Q)  satisfying  Assumptions 
4. 1-4.4.  Let  C(s)  A  Ne(s)/De(s).  It  follows  from  Corollary  4.1  that  we  need 
to  show 

A(s,q)  A  Dc(s)D(s ,q)  +  Ne(s)N(a,q)  (J.l) 

is  Hurwitz  invariant  over  Q  and 

_ |Nf(»N(JM,q)  |  | r(ja>)  | 

|De(j«)D(jW‘q)  +  Hc<J«)N(j«,q)  iP6(Jw,q)  |  ~ 


for  all  q  e  Q  and  w  >  0.  By  Assumption  4.4,  we  have 

»l 

in 


r  |r(J«) 

sup[ 


q  €  Q 


;  «  >  o| 


<  l. 


Consequently,  there  exists  a  S  >  0  such  that 

T.lfflSjU  <  1/(1  ♦  «)■ 

for  all  q  €  Q  and  u  >  0.  Hence,  we  only  need  to  show  that 

iN.QOHQ.q)  |  .  z  ,  s 

|De(J»)D(j«,q)  +  Ne(jw)N(Jw,q)  |  - 


(J.2) 


for  all  q  €  Q  and  u  >  0. 

To  construct  a  compensator  C(s)  satisfying  (J.l)  and  (J.2),  begin  by 
choosing 

k 

Nc(s)  A  £  V1 
i-0 

to  be  any  Hurwitz  polynomial  with  \  >  0  when  N(s,q)  is  PHI,  or  X*  <  0  when 
N(s,q)  is  negative  Hurwitz  invariant  (NHI)  and  of  degree 
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k  >  max  (d'(q)  -  d(q) :  q  €  Q)  -  1. 

Note  that  if  max  (d'(q)  -  d(q) :  q  €  Q|  -  0,  we  simply  take  k  >  0. 

Recall  that  D(s,q)  is  a  standard  uncertain  polynomial  and  N(s,q)  Is  a 
standard  PHI  or  NHI  polynomial.  It  follows  that  A^s.q)  A  Nc(s)N(s,q)  Is  also 
a  standard  PHI  polynomial.  Take 
f(s,q)  -  A0(s , q)  and 
f ' (s,q)  -  D(s , q) . 

It  is  easy  to  verify  that  all  hypotheses  of  Lemma  4.2  are  satisfied.  Hence, 
we  can  select  a  sufficiently  small  cQ  >  0  such  that 
AjCs.q)  A  A0(s,q)  +  c0D(s,q) 
is  a  standard  PHI  polynomial  and 

<  1  +  S/2 

for  all  q  €  Q  and  u  >  0.  The  proof  la  continued  by  repeated  application  of 
Lemma  4.2;  i . e . ,  take 
f (s  ,q)  -  A^s.q) 
f'(s,q)  -  sD(s,q) 

and  again,  the  requirements  of  the  lemma  are  satisfied.  Hence,  we  can  choose 

cx  >  0  small  enough  that 

A2(s,q)  A  A0(s,q)  +  (c0  +  Cj,s)D(a,q) 

is  a  standard  PHI  polynomial  and 

|A1<Jw,q)t  1  +  S/2  +  */4 
|Aj<J«.q)  |  ~  1  +  S/2 

for  all  q  €  Q  and  w  >  0.  Continuing  in  this  manner,  we  now  apply  Lemma  4.2 
with  f(s,q)  -  Aj(s,q)  and  f'(s,q)  -  s3D(s,q),  and  can  generate  c2  >  0  such 
that 
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lA^jw.q)  |  |A2(jw,q)  |  ...  |Aj(j»,q)| 

<  1  +  6/2  +  6/4  +  ...  +  6/2*  <1+6 
for  all  q  €  Q  and  <■>  >  0. 

A  similar  proof  applies  to  an  uncertain  system  M(P0(s,q) , r(s) ,Q) . 
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